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ABSTRACT 

We report the results of a search for 12.2-GHz methanol maser emission, targeted 
towards 113 known 6.7-GHz methanol masers associated with 1.2-mm dust continuum 
emission. Observations were carried out with the Australia Telescope National Facility 
(ATNF) Parkes 64-m radio telescope in the period 2008 June 20 - 25. We detect 68 
12.2-GHz methanol masers with flux densities in excess of our 5-cr detection limit of 
0.55 Jy, 30 of which are new discoveries. This equates to a detection rate of 60 per 
cent, similar to previous searches of comparable sensitivity. 

We have made a statistical investigation of the properties of the 1.2-mm dust 
clumps with and without associated 6.7-GHz methanol maser and find that 6.7-GHz 
methanol masers are associated with 1.2-mm dust clumps with high flux densities, 
masses and radii. We additionally find that 6.7-GHz methanol masers with higher peak 
luminosities are associated with less dense 1.2-mm dust clumps than those 6.7-GHz 
methanol masers with lower luminosities. We suggest that this indicates that more 
luminous 6.7-GHz methanol masers are generally associated with a later evolutionary 
phase of massive star formation than less luminous 6.7-GHz methanol maser sources. 

Analysis of the 6.7-GHz associated 1.2-mm dust clumps with and without asso- 
ciated 12.2-GHz methanol maser emission shows that clumps associated with both 
class II methanol maser transitions are less dense than those with no associated 12.2- 
GHz methanol maser emission. Furthermore, 12.2-GHz methanol masers are prefer- 
entially detected towards 6.7-GHz methanol masers with associated OH masers, sug- 
gesting that 12.2-GHz methanol masers are associated with a later evolutionary phase 
of massive star formation. 

We have compared the colours of the GLIMPSE point sources associated with 
the maser sources in the following two subgroups: 6.7-GHz methanol masers with and 
without associated 12.2-GHz methanol masers; and 6.7-GHz methanol masers with 
high and and those with low peak luminosities. There is little difference in the nature 
of the associated GLIMPSE point sources in any of these subgroups, and we propose 
that the masers themselves are probably much more sensitive than mid-infrared data 
to evolutionary changes in the massive star formation regions that they are associated 
with. 

We present an evolutionary sequence for masers in high-mass star formation re- 
gions, placing quantitative estimates on the relative lifetimes for the first time. 
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1 INTRODUCTION 

Interstellar masers are one of the most readily observed 
signposts of star formation, particularly emission from 
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the hydroxyl (OH), water (H 2 0) and methanol (CH 3 OH) 
molecules as they are relatively common and intense. Be- 
cause these masers arise at radio frequencies, they are not 
affected by the dense optical-obscuring gas and dust present 
at the early stages of massive star formation thereby allow- 
ing us to probe these stars at the earliest stages. The masers 
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they are known to typ ically be co-spatial (|Menten et al.l 



are powerful probes of the kinematics of the star formation 
regions, but we are only recently beginning to be able to use 
the masers to investigate aspects such as the physical condi- 
tions through comp arison of observat i onal data and mase r 
pumping theories (ICragg et al.1 l200ll ; ISutton et ail l200lf ). 
The different maser species favour different physical condi- 
tions and are thought to trace different evoluti onary phases 
of massive star formation (jEllingsen et al.ll2007l ). since many 
sources show emission from multiple maser transitions or 
species there must be significant overlap for the evolution- 
ary phase traced by the most common types of masers. 

Interstellar masers have been detected towards more 
than a thousand star formation regions in our Galaxy from 
transitions of methanol, hydroxyl and water. Transitions of 
class II methanol masers have some advantages over water 
and OH m asers as they exclusively trace sites o f massive star 
formation (|Minier et al.l 120031 ; IXu et alj I200ST ) while water 
and OH are known to also be associated with other astro- 
physical objects such as evolved stars and supernova rem- 
nants. Two of the strongest maser species found towards star 
formation regions are transitions of methanol at 6.7- and 
12.2-GHz, both of which are known to tr ace an early evo - 
lutionary stage of massive star formation l|Ellingsenl 2006). 
Complementary observations of these two strong methanol 
maser transitions are especially useful in probing the physi- 
cal conditions of the environments in which they arise, as 

;yp i< 

ll992l ; lNorris et alj|l993h . 

To date there have been several searches for 12.2-GHz 
methanol maser emission towards 6.7 - GHz methanol maser 
emiss i on (e.g. | Caswell et al.l Il995bl ; iBIaszkiewicz fe Kui 
l2004MCasweH et al.1 (|l995bh targeted their search for 12.2- 
GHz methanol masers towards 238 6.7-GHz methanol 
masers that had been detected towards sites of OH 
masers and achieved a d etection rate of 55 per cent. 
Blaszkicwicz & Kul I 20041 ) directed their observations at 
12.2-GHz towards 6.7-GHz methanol masers that were ob- 
served towards a sample of IRAS (Infrared Astronomy 
Satellite) selected sources l|Szvmczak et alJl200Ch as well as 
6.7 GHz methanol masers that were detected in a blind 
search t hat had relatively poo r positional accuracy and sen- 
sitivity l|Szvmczak et alj|2002h . 

Searches such a s those by ICaswell et al.l l| 1995b! ) and 
IBIaszkiewicz fc Kud l|2004h have shown that 12.2-GHz 
methanol maser emission is only rarely brighter than the 
associated 6.7-GHz methanol maser emission. This coupled 
with the fact that there have been no serendipitous detec- 
tions of 12.2-GHz masers without a 6.7-GHz counterpart 
mean that it is unlikely that an unbiased Galactic survey 
would uncover many, if any, more 12.2-GHz methanol masers 
than would be yielded from a search targeted towards 6.7- 
GHz methanol maser sources. 

Recently, a small subset of the s trong 6.7-GHz m ethanol 
masers from the sources observed bv lHill et al.l l|2005h at 1.2- 
mm, were searched for the presence of 12.2-GHz methanol 
masers with the Uni versity of Ta s mania 26 m M ount Pleas- 
ant radio telescope (|Lewisl |200"7T I. iLewisI l|2007h detected 17 
12.2-GHz methanol masers tow ards the 27 sites searched 
with a detection limit of ~4 Jy. iHill et all (120051 ) observed 
131 regions that were suspected of undergoing massive star 
formation, using the presence of previously identified 6.7- 
GHz methanol masers and/or an ultracompact Hn region 



to select these regions. A total of 404 1.2-mm dust clumps 
were identified, with 113 of these associated with 6.7-GHz 
methanol masers, 35 of which are as sociated with radio con- 
tinuum emission (|Walsh et al. 19981). 

Statistical analysis by Lewis! (|2007) of the properties 
of 1.2-mm dust clumps associated with 6.7-GHz methanol 
maser emission with and without 12.2-GHz methanol masers 
showed that those 1.2-mm dust clumps with both 6.7-GHz 
and 12.2-GHz were less massive and less dense than the 
1.2-mm clumps that were associated with only 6.7-GHz 
methanol maser emission. Comparison between the 6.7-GHz 
methanol masers with detected 12.2-GHz with those with 
no detectable 12.2-GHz maser emission, showed that 12.2- 
GHz methanol masers were preferentially found at those 
6.7-GHz methanol maser sites with associated OH emission 
(to a better than 95 per cent confidence). As OH masers 
are e xpected to trace a generally later stage of star forma- 
tion (|Forster fc Caswelllll989l ; ICaswelll Il997l ) than 6 . 7-GHz 
methanol masers. As the sample size of the ILewisI (|2007D 
observations was small, was subject to a 6.7-GHz flux den- 
sity bias and had relatively poor sensitivity, the aim of the 
current work was to confirm these results by testing a large 
sample with a much lower 12.2-GHz detection limit. 

Here we present a targeted search for 12.2-GHz 
methanol mase rs, wit h a gr eatly increased sample size and 
sensitivity (c.f. iLewid (|2007l )). towards 113 known 6.7-GHz 
methanol masers which have been previously targ eted for 
1.2-mm dust continuum emission (Hill ct al. 2005). Statis- 
tical analysis of the 1.2-mm dust clumps devoid of either 
methanol maser transition and those associated with 12.2- 
GHz and/or 6.7-GHz methanol maser sources has been car- 
ried out, with a particular emphasis on the investigation of 
the relative evolutionary stage each class of source is tracing. 



2 OBSERVATIONS & DATA REDUCTION 

We h ave searched 113 6.7-GHz methanol masers, towards 
which IHill et al.1 (2005) observed 1.2-mm dust clump con- 
tinuum emission, for the presence of associated 12.2-GHz 
methanol maser emission. Observations were carried out 
with the Australia Telescope National Facility (ATNF) 
Parkes 64-m radio telescope during 2008 June 20 - 25. 
Follow-up observations of two sources were made on 2008 
December 7 using an identical set-up to the observations car- 
ried out in 2008 June. The observations were made with the 
Ku-band receiver which detected two orthogonal linear po- 
larizations and had typical system temperatures of 205 and 
225 Jy for the respective polarizations throughout the obser- 
vations. The Parkes multibeam correlator was configured to 
record 8192 channels over 16-MHz for each of the recorded 
linearly polarized signals. This configuration yielded a use- 
able velocity coverage of ~290 km s _1 and a spectral reso- 
lution of 0.08 km s _1 , after Hanning smoothing. The Parkes 
radio telescope has rms pointing errors of ~10 arcsec and at 
12.2-GHz the telescope has a half power beam width of 1.9 
arcmin. 

All sources were observed at a fixed frequency of 12178- 
MHz (i.e. with no Doppler tracking), which alleviated the 
requirement for a unique reference observation to be made 
for each of the source positions. Instead, all the spectra col- 
lected each day were combined to form a sensitive reference 
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bandpass by taking the median value of each channel. The 
noise contribution to the quotient spectrum from such a ref- 
erence is negligible and means that the theoretical radiome- 
ter noise is achieved for each observation. 

As the bandpass is not constant throughout an observ- 
ing session, observing in this mode introduced ripples into 
the quotient spectra. This was overcome by subtracting a 
running median over 100 channels from the baseline. As the 
maser lines are much narrower than 100 channels this has lit- 
tle effect on emission spectra but may have an adverse effect 
on absorption features which are of a comparable width. Our 
ability to detect absorption is therefore reduced and in most 
cases where we do detect absorption it appears narrower and 
weaker than it would have been if we had employed a more 
traditional, but less efficient, observing method. In order 
to determine the absorption characteristics of the sources, 
spectra where absorption features were detected were box- 
car smoothed over 20 channels, decreasing our spectra reso- 
lution to 0.26 km s _1 and resulted in a typical rms noise of 
-0.04 Jy. 

Data were reduced using ASAP (ATNF Spectral Anal- 
ysis Package). Alignment of velocity channels was carried 
out during processing. Absolute flux density calibration was 
achieved by observing PKS B1934-638 each day whi ch has 
an as sumed flux density of 1.825 Jy at 12178-MHz (|Saultl 
2003). Each target was observed for 10 minutes, and once 
the polarizations were averaged yielded a typical rms noise of 
0.11 Jy corresponding to a 5-a detection limit of 0.55 Jy. For 
weak sources (of comparable strength to our detection limit), 
individual polarizations were inspected as a further check on 
the reliability of the detection. The adopted rest fre quency 
was 12.178597 GHz (jMiiller. Menten fc Maderll2004h . 



therefore our observations targeted a total of 106 1.2-mm 
dust clumps. Spectra of the detected 12.2-GHz methanol 
masers are shown in Fig. Q] and are presented in order of 
Galactic longitude except for where vertical alignment of 
nearby sources was necessary to highlight where features of 
a source were also detected at nearby positions. Comments 
on some individual sources can be found in Section [3TT| The 
absorption features that we detect are presented in Table [2] 



3 RESULTS 



Our search, for 12.2-GHz methanol masers, targeted 
113 6.7-GHz methanol maser sites and resulted in 
the detection of 68 12.2-GHz methanol masers, equat- 



of the literature (e.g. 


Pestalozzi. Minier & Boothl 12005; 


Blaszkicwicz & Kus 2004 


;ICaswell et al.lll995bh reveals that 



30 of these are new detections. A number of the 12.2- 
GHz methanol masers that have been previously ob- 
served, particularly the stronger ones, were di s covered 
prior to iBIaszkiewicz fc Kusl i2004): ICaswell et all 111995b ) 



by others (e.g. | Gavlard. MacLeod fc van der Waltl 11994 ; 
ICaswell et al.lll993h . However, direct comparison of the ma- 
jority of these sources are complicated by large positional 
uncertainties, pointing offsets and the limited velocity ranges 
covered in the observations. 

Table [T] summarizes in formation on the targeted 
6. 7- GHz methanol masers llCaswelll 120091; Minier et al. 
20031; IWalsh et all Il998l, Il997l ; ICaswell et al.l Il995, 



Szvmczak et al.l 20021 ) . associated 12.2-GHz methanol 



masers that we detect, along with the 1.2-mm dust clump 
(|Hill et al.ll2005l ) that each methanol maser site is associ- 
ated with. For the majority of sources, precise positions for 
the 6.7-GHz methanol maser sources have been determined 
with the Australia Telescope Compact Array to an accu- 
racy better than 1 arcsec. There are five instances where two 
methanol masers are associated with the same dust clump, 
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Table 1. Targeted 6.7-GHz methanol masers, detected 12.2-GHz methanol masers and associated 1.2-mm dust clump that they are 
associated with. Column 1 gives the Galactic coordinates for each methanol maser source, derived from the position of the 6.7-GHz 
methanol maser, columns 2 and 3 give the 6.7-GHz methanol maser right ascension and declinati on. References for each 6.7-GHz met hanol 
maser position are rep resented by the sup ersc ript number at the end of the source name where 1 : ICaswelll (I2009T) , 2: lMinier et alJd2003ft . 3: 
IWalsh et aj]<199Sl) . 4: IWalsh et all dl997l ). 5: ICasweII et al.l dl995al) , 6: ISzvmczak et alj ll2002h . a: James Caswell, private communication; 
see also note in comments on sources for G 034.246+0.134 and G 035.025+0.350. References for previously detected 12.2-GHz methanol 
masers are listed after the 6.7-GHz met han ol maser references a n d are separated by a comma and contained within brackets to avoid 
confusion, here 7: ICaswell et al ] dl995bh . 8: iBlaszkiewicz fc Kusl j2004h . *: new detection. Colu mns 4-10 list the 6.7-GHz peak flux 
density (Jy), peak ve locity (km s — 1 ), vel ocity range (km s — 1 ) (nb. where the source r eference is lWalsh et al. I ll 19981) the velocity range 
has been taken from IWalsh et al ] l| 19971 ) where available and Pcstalozzi ct al . (2005) if not), 12.2-GHz peak flux density (Jy), peak 
velocity (km s _1 ), velocity range (km s — 1 ) and integrated flux density (Jy km s — 1 ). The 5-<r detection limit and velocity range covered 
at 12.2-GHz is presented for sources where we detect no 12.2-GHz emission, in these cases the 5-<r detection limit is presented in column 
7 and is preceded by a < , and the velocity ranges observed are listed in column 9. @ indicates a source that was previously observed 
bv lCaswell et al] l|l995bl) but we do n o t dete ct. Column 11 indicates if the methanol maser sites have been searched for the presence 
of an associated OH maser bv lCaswelll d 19981) . — indicates that no suitable data were found, n indicates that the site has been searched 
but no source w as detected and y indicates that there is an OH maser within 2 arcsec of the 6.7-GHz methanol maser. The dust clump 
dHill et al.ll2005l) that the methanol masers are coincident with are listed in column 12, a * after the source name indicates sources that 
were not included in the statistical analysis, for details see Section l3.2l 
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11 


26 


59 


-51 


41 


56.7 


190 


-78.2 


-87,-77 


76 


-78.7 


-83.1,-77.7 


58.8 


y 


G 331.28-0.19* 


G332.653-0.621 1 


16 


19 


43 


51 


-51 


03 


39.9 


7.1 


-50.6 


-52,-49 


<0.54 




-130,158 




ii 


G 332.648-0.606 


G332.726-0.621 1 '(*' 


16 


20 


03 


00 


-51 


00 


32.5 


2.7 


-49.5 


-57,-44 


0.5 


-47.3 


-47.4,-47.2 


0.1 


y 


G 332.73-0.62 


G000.212-0.001 1 '(*) 


17 


46 


07 


63 


-28 


15 


20.9 


3.5 


49.2 


41,50 


0.8 


49.5 


49.1,49.7 


0.5 


n 


G 0.21-0.00 


GOOO.SIS-O^OI 1 '**) 


17 


47 


09 


13 


-28 


46 


15.7 


41.2 


18 


14,27 


2.4 


18.3 


17.8,18.5 


1.0 


n 


G 0.32-0.20 


G 000.316-0. 201 1 


17 


47 


09 


33 


-28 


46 


16.0 


1.3 


21 


20,22 


<0.55 




-120,168 




ii 


G 0.32-0.20 
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Table 1. - continued 



Methanol maser 

(l,b) 
(degrees) 


RA 

(J2000) 
(h m s) 


Dec 
(J2000) 

(° ' ") 


£>6.7 

(Jy) 


Vpg.7 


Vr 6 . 7 


Si2.2 

(Jy) 


Vpi2.2 


Vri2.2 


Il2.2 


OH 


Dust clump 

(l,b) 

(degrees) 


G 000.496+0. lSS 1 ^*) 


17 46 03.96 


-28 24 52.8 


10 


0.8 


-12,2 


18 


1.1 


-9.5,1.5 


15.6 


y 


GO.49+0.19 


G000.546-0.852 1 >< 7 > 


17 50 14.35 


-28 54 31.1 


50 


13.8 


8,20 


3.9 


18.5 


16.1,19 


2.1 


y 


G 0.55-0.85 


G000.836+0.184 1 


17 46 52.86 


-28 07 34.8 


8.1 


3.5 


2,4 


<0.52 




-115,173 




n 


G 0.83+0.18 


G 001. 147-0. 124 3 


17 48 48.467 


-28 01 11.81 


4 


-20.8 


-21,-15 


<0.51 




-115,173 




n 


G 1.14-0.12 


G002.536+0.198 1 'M 


17 50 46.47 


-26 39 45.3 


40 


3.2 


2,20 


3.7 


7.0 


3.4,7.2 


2.6 


n 


G 2.54+0.20 


G005.900-0.430 1 


18 00 40.86 


-24 04 20.8 


5.3 


10 


4,11 


<0.55 




-113,175 




n 


G 5.90-0.42 


G 006.610-0. 082 1 >M 


18 00 54.03 


-23 17 03.1 


10.8 


0.7 


0,1 


10.7 


0.9 


0.5,1.4 


4.7 


n 


G 6.60-0.08* 


G 008.139+0. 226 1 


18 03 00.75 


-21 48 09.9 


3.5 


20 


19,21 


<0.55 




-112,176 




n 


G8. 13+0. 22 


Goos.eeg-o.sse 1 ^*) 


18 06 18.99 


-21 37 32.2 


8.5 


39.3 


39,40 


0.3 


39.5 


39.4,39.6 


0.05 


y 


G 8.68-0.36 


Goos.ess-o.ses 1 '^ 


18 06 23.49 


-21 37 10.2 


70 


42.9 


40,46 


8.4 


44.2 


41.0,45.5 


11.6 


y 


G 8.686-0.366 


G 009.619+0. 193^ 


18 06 14.92 


-20 31 44.3 


72 


5.5 


5,7 


0.4 


5.8 


5.7,5.9 


0.1 


y 


G 9.62+0.19 


G 009.621+0. WG 1 '^ 


18 06 14.67 


-20 31 32.4 


5000 


1.3 


-4,9 


401 


1.4 


-3.0,5.8 


202.8 


y 


G9. 62+0. 19 


G009.986-0.028 1 '(*> 


18 07 50.12 


-20 18 56.5 


28 


47.1 


40,52 


1.0 


47.1 


44.3,47.3 


0.5 


n 


G 9.99-0.03 


G 010.287-0. 125 1 


18 08 49.36 


-20 05 59.0 


27 


5.0 


4,6 


<0.52 




-111,178 




n 


G 10.284-0.126 


G 010.299-0. 146 1 


18 08 55.54 


-20 05 57.5 


3.5 


20 


19,21 


<0.54 




-111,178 




n 


G 10.29-0.14 


G 010. 323-0. 160 1 ''*) 


18 09 01.46 


-20 05 07.8 


126 


10 


4,14 


2.3 


10.3 


5.1,12.6 


1.8 


n 


G 10.32-0.15 


G 010.444-0. OlS 1 '' 7 ) 


18 08 44.88 


-19 54 38.3 


14.8 


73.2 


68,79 


6.5 


72.0 


68.0,73.7 


10.1 


y 


G 10.44-0.01 


G 010.473+0. 027 1 '( 7 > 


18 08 38.20 


-19 51 50.1 


120 


75.0 


58,77 


12.4 


75.1 


73.7,76.7 


9.1 


y 


G 10.47+0.02 


G 010.480+0. 033 1 


18 08 37.88 


-19 51 16.0 


9.7 


65.0 


58,66 


<0.55 




-111,178 




y 


G 10.47+0.02 


G010.627-0.384 1 '( 7 > 


18 10 29.22 


-19 55 41.1 


3.1 


4.6 


-6,7 


1.4 


4.6 


4.0,6.6 


0.8 


n 


G 10.62-0.38 


G010.629-0.333 1 


18 10 17.98 


-19 54 04.8 


4.2 


-7.5 


-13,1 


<0.54 




-111,178 




n 


G 10.62-0.33 


G 011.497-1. 485 1 -M 


18 16 22.13 


-19 41 27.1 


167 


6.7 


4,17 


29 


9.0 


6.0,17.1 


15.8 


n 


G 11.49-1.48 


G011.904-0.141 1 '( 7 ) 


18 12 11.44 


-18 41 28.6 


56 


42.8 


40,45 


17 


43.0 


42.1,44.0 


8.2 


y 


G 11.903-0.140 


G 011. 936-0. 150 1 


18 12 17.29 


-18 40 02.6 


1.9 


48.5 


47,50 


<0.76 




-111,178 




n 


G 11.93-0.14 


G 011. 936-0. eie 1 '' 7 ) 


18 14 00.89 


-18 53 26.6 


41 


32.2 


30,44 


4.3 


39.8 


31.5,42.4 


4.4 


n 


G 11.93-0.61 


G011.992-0.272 3 


18 12 51.197 


-18 40 39.7 


0.6 


59.5 


59,61 


<0.75 




-110,179 




n 


G 11.99-0.27 


G 012.025-0. 025 1 '( 7 > 


18 12 01.86 


-18 31 55.7 


85 


107.7 


105,112 


14 


108.6 


107.8,109.4 


8.8 


n 


G 12.02-0.03 


G 012.202-0. 120 1 


18 12 42.93 


-18 25 11.8 


1.7 


26.4 


26,27 


<0.75 




-111,178 




n 


G 12.18-0.12 


G 012.209-0. 102 1 


18 12 39.92 


-18 24 17.9 


9.2 


19.8 


16,22 


<0.75 




-111,178 




n 


G 12.20-0.09 


G 012.681-0. 182 1 >( 7 ) 


18 13 54.75 


-18 01 46.6 


450 


57.6 


50,61 


8 


57.1 


55.6,60.3 


7.6 


y 


G 12.68-0.18 


G012.889+0.489 1 '( 7 > 


18 11 51.39 


-17 31 29.6 


27 


39.3 


28,42 


19 


39.4 


38.9,40.2 


8.7 


y 


G 12.88+0.48 


G012.909-0.260 1 '( 7 > 


18 14 39.53 


-17 52 00.0 


300 


39.9 


35,47 


16 


39.4 


38.5,40.6 


14.0 


y 


G 12.90-0.26 


G014.60+0.01 3 'M 


18 17 01.139 


-16 14 38.88 


1.2 


23.2 


22,28 


0.6 


24.7 


24.6,24.7 


0.1 


n 


G 14.60+0.01 


G015.034-0.677 1 '( 7 > 


18 20 24.78 


-16 11 34.6 


28 


21.2 


20,24 


21 


23.4 


21.2,23.7 


12.1 


y 


G 15.03-0.67 


G016.58-0.05 3 '( 7 > 


18 21 09.132 


-14 31 48.69 


22 


59.3 


52,69 


2.2 


65.9 


59.2,68.7 


2.2 


y 


G 16.58-0.05 


G016.86-2.15 3 


18 29 24.411 


-15 16 04.12 


24 


14.8 


14,19 


<0.53 




-108,180 




y 


G 16.86-2.15 


G 019.609-0.234" 


18 27 37.99 


-11 56 37.6 


0.4 


40.1 


36,43 


<0.53 




-107,181 




y 


G 19.607-0.234 


G019.612-0.134 3 '( 7 > 


18 27 16.524 


-11 53 38.49 


18 


56.2 


50,61 


1.3 


50.5 


50.2,55.1 


0.8 


- 


G 19.61-0.1 


G021.880+0.014 3 


18 31 01.749 


-09 49 01.13 


5.5 


20.6 


17,22 


<0.52 




-106,182 




- 


G 21.87+0.01 


G 022.356+0. 066 3 'M 


18 31 44.128 


-09 22 12.40 


10.0 


80.0 


77,88 


1.0 


81.4 


81.4,81.5 


0.1 


- 


G 22.35+0.06 


G023.257-0.241 3 >M 


18 34 31.259 


-08 42 46.84 


6.0 


64.0 


63,66 


0.9 


64.9 


63.8,65.9 


0.7 


- 


G 23.25-0.24 


G 023.437-0. 184 1 '( 7 > 


18 34 39.25 


-08 31 38.5 


45 


103.0 


101,108 


3.3 


104.0 


103.6,105.3 


1.4 


y 


G 23.43-0.18 


G023.440-0.182 1 '( 7 > 


18 34 39.18 


-08 31 25.3 


25 


96.6 


94,100 


9 


97.5 


96.1,107.3 


16.0 


- 


G 23.43-0.18 


G024.790+0.083 3 


18 36 12.571 


-07 12 11.31 


63 


113.3 


106,115 


<0.52 




-104,184 




y 


G 24.78+0.08 


G 024.850+0. 087 3 -(*> 


18 36 18.398 


-07 08 52.09 


24 


109.7 


109,115 


1.8 


110.1 


107.7,110.7 


1.2 


- 


G 24.84+0.08 


G 025.650+1. 050 3 


18 34 20.913 


-05 59 40.44 


90 


41.7 


38,43 


<0.51 




-104,184 




- 


G 25.65+1.04 


G 025.710+0. 044 3 >( g ) 


18 38 03.148 


-06 24 14.80 


355 


95.5 


89,101 


144 


95.6 


89.3,100.4 


107.7 


- 


G 25.70+0.04 


G025.826-0.178 3 'M 


18 39 03.631 


-06 24 09.64 


36 


91.5 


89,94 


25 


90.9 


90.0,92.1 


14.7 


- 


G 25.83-0.18 


\j 025. 14b— U.OUo A ' 


lo 4z 4Z.D9 


C\A 1 K Id 

—04 15 oD.O 


61 


101.2 


99,105 


20 


101.3 


no o i r\ a 1 
y9.y, 104.1 


11.9 




(~* OO 1 A f\ C\f\ 

Lt 28.14—0.00 


G028.201-0.049 1 


18 42 58.08 


-04 13 56.2 


3.5 


98.9 


94,100 


<0.52 




-103,185 




y 


G 28.20-0.04 


G028.282-0.359 3 '(*> 


18 44 13.270 


-04 18 02.99 


1.1 


40.4 




0.7 


41.0 


40.8,41.0 


0.1 




G 28.28-0.35 


G028.305-0.387 3 'M 


18 44 21.988 


-04 17 38.42 


66 


81.0 


79,93 


8 


81.9 


80.8,83.0 


7.9 




G 28.31-0.38 


G029.865-0.043 3 '( 7 > 


18 45 59.588 


-02 44 57.60 


12 


101.2 


95,105 


21 


100.4 


99.9,102.1 


10.1 




G 29.86-0.04 


G029.956-0.016 3 '( 7 > 


18 46 03.763 


-02 39 19.69 


158 


95.7 


95,102 


51 


96.6 


91.2,98.5 


45.0 




G29.96-0.02B 


G 029.979-0. 047 3 '( 7 > 


18 46 12.974 


-02 38 58.05 


14 


98.2 


97,105 


6.6 


98.3 


96.5,103.9 


4.9 




G 29.978-0.050 


G030.591-0.042 3 


18 47 18.886 


-02 06 07.00 


0.6 


43.2 


36,49 


<0.55 




-102,186 






G 30.59-0.04 


G030.704-0.068 3 


18 47 36.801 


-02 00 48.99 


66 


88.0 


85,90 


<0.55 




-102,186 






G 30.705-0.065 


G 030.760-0. 052 3 '( 7 > 


18 47 39.729 


-01 57 21.92 


75 


91.8 


88,95 


1.7 


91.3 


90.7,91.8 


0.9 




G 30.769-0.048 


G030.818-0.057 3 '( 7 > 


18 47 46.984 


-01 54 19.61 


10 


101.2 


91,110 


0.5 


108.1 


108.0,108.1 


0.1 




G 30.81-0.05 


G030.898+0.162 3 'M 


18 47 09.133 


-01 44 10.29 


40 


101.8 




14 


101.7 


99.3,102.1 


13.5 




G 30.89+0.16 
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Table 1. - continued 



lV/Tci+Viiinnl YnncoT* 


RA 


Dec 


§6.7 


Vp6.7 


Vre i 
vr 6.7 


S12.2 


Vpi2.2 


Vri2.2 


I12.2 


OH 




(l,b) 


(J2000) 


(J2000) 


(Jy) 






(Jy) 










(l,b) 


(degrees) 


(h m s) 


(° ' ") 


















(degrees) 


G031.061+0.094 3 


18 47 41.342 


-01 37 21.31 


16 


16.2 


15,17 


<0.5 




-102,186 






G 31.06+0.09 


G 031. 282+0. 062 3 >< 7 ) 


18 48 12.385 


-01 26 22.60 


83 


110.0 


102,113 


15 


110.4 


105.1,113.0 


25.8 




G 31.28+0.06 


G031.411+0.307 3 '( 7 > 


18 47 34.314 


-01 12 47.14 


9 


103.4 


90,108 


1.5 


95.7 


94.1,95.8 


0.6 


y 


G 31.41+0.30 


G033.13-0.09 5 '( 7 > 


18 52 07.2 


+00 08 06 


12.4 


73 


71,81 


0.35 


73.3 


72.9,73.3 


0.1 




G 33.13-0.09 


G034.246+0.134 a >( 7 ) 


18 53 21.46 


+01 13 54 


10.6 


55.4 


55,63 


1.3 


61.2 


61.0,61.8 


0.8 




G 34.24+0.13 


G 035.025+0.350"' W 


18 54 00.66 


+02 01 18 


36 


44.3 


40,47 


0.7 


42.2 


42.0,43.4 


0.2 




G 35.02+0.35 


G 037.47-0. ll 6 'M 


19 00 06.7 


+03 59 27 


12.5 


58.4 


54,64 


2.0 


63.0 


58.6,63.1 


1.2 




G 37.475-0.106 


G049.470-0.371 1 'M 


19 23 37.90 


+14 29 59.4 


9.2 


64.0 


55,76 


2.0 


57.8 


57.0,64.0 


1.2 




G 49.472-0.366 


G049.489-0.369 1 '( 7 > 


19 23 39.83 


+14 31 05.0 


26 


56.1 


55,61 


4.8 


56.2 


55.6,57.3 


2.1 




G 49.49-0.37 
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3.1 Comments on individual sources 

G 192. 600-0. 048 : The weak 12.2-GHz methanol maser 
emission towards this 6. 7- GHz methanol maser is one of the 
two weakest of our detections. The emission has a peak flux 
density of 0.47 Jy, which corresponds to a 4.7-a detection. 
Furthermore, the velocity of the 12.2-GHz emission is 3.6 
km s~\ close to the peak of the 6.7-GHz methanol maser at 
5 km s _1 , well within the full velocity range of the associated 
6.7-GHz methanol maser. 

G 284.352-0.419 : ICaswell etahl i|l995bl ) detected a 
0.7 Jy 12.2-GHz methanol maser towards this 6.7-GHz 
methanol maser in 1992. We detect no hint of emission at 
12.2-GHz and expect that this is due to source variability. 
Due to this we have excluded this source from our statistical 
analysis. 

G 290.374 + 1.661 : Observations in 2008 June detected 
marginal 12.2-GHz emission toward this 6.7-GHz methanol 
maser outside the velocity range of the 6.7-GHz maser emis- 
sion at a peak velocity of -35.1 km s~ of 0.6 Jy. Follow-up 
observations carried out in 2008 December uncovered no de- 
tectable 12.2-GHz emission. We have therefore listed this 
source as having no associated 12.2-GHz emission but have 
removed this source from our statistical analysis to avoid 
the possibility that we have misclassified this source. 

G 291.579-0.431 : During the 2008 June observations 
we detected weak 12.2-GHz methanol maser emission to- 
wards this source. The velocity of the detected emission was 
within the velocity range of the 6.7-GHz methanol maser 
emission, however, since the detection was less than 5-cr the 
source was re-observed in 2008 December and resulted in 
no detectable emission. As for G 290.374+1.661, we have re- 
moved this source and the associated 1.2- mm dust clump 
from our statistical analysis. 

G 305.208+0.206 : Observations towards 6.7-GHz 
methanol maser G 305.208+0.206 as well as near-by 
G 305.202+0.208, offset by 22 arcsec, were carried out in ad- 
jacent observations. The relative amplitudes of the detected 
emission at these two positions suggest that all of the emis- 
sion is associated with G 305.208+0.206; however, the peak 
of the 12.2-GHz emission is at the p eak velocity of nearby 
G 305.202+0.208 and observations bv lCaswell et alJ i|l995bh 
also indicate that there is 12.2-GHz emission associated with 
each 6.7-GHz methanol maser source. We suggest that either 
pointing errors or unfavourable weather conditions have re- 
sulted in the relative amplitudes of the emission towards the 
two sources being affected. Both of these methanol maser 
sites are located within the one dust clump. 

G 305.366+0.184, G 000.496+0.188, G 029.865-0.043 : 
These sources all have stronger 12.2-GHz methanol maser 
emission than the values for their 6.7-GHz counterparts 
as shown in Table [1] As masers are intrinsically variable, 
measurements of the 6.7-GHz methanol maser flux density 
would need to be taken at the same epoch as the 12.2-GHz 
methanol masers to be certain that this was not a result of 
source variability. 

G 034.246+0.134 and G 035.025+0.350 : The positions 
for these sources were determined with the ATCA with a 
large declination uncertainty of ±15 arcsec (James Caswell, 
private communication) owing to the highly elongated beam 
for sources near declination 0°. 

G 049.489-0.369 : We observed this source towards a 



position offset from the true position by ~66 arcsec (58 
percent of the HPBW). We therefore expect that the flux 
density of this source is underestimated by approximately a 
factor of two. 

3.2 Statistical analysis 

In order to determine if there are any differences between 
1.2-mm dust clumps without associated 6.7-GHz methanol 
masers, those associated with 6.7-GHz methanol masers only 
and those associated with both 6.7- and 12.2-GHz methanol 
maser emission, we c arried out a statistic al analysis following 
the method used by iBreen et all (|2007i ) to investigate wa- 
ter maser presence towards 13 CO clumps and 1.2-mm dust 
clumps. Additionally, the possibility of a relationship be- 
tween dust clump properties and OH maser presence as well 
as peak 6.7-GHz methanol maser luminosity values and the 
presence of detectable radio continuum were investigated. 
Due to the difficulty in determining meaningful uncertain- 
ties for many of the tested source characteristics, all values 
have been equally weighted and assumed to be error free in 
all subsequent statistical analysis. 

Of the 40 4 1.2-mm dust clump sources listed by 
iHill et al.li|2005h . 26 were excluded from our statistical anal- 
ysis as they were either missing a value for one or more 
of the dust c lump properties to be tested, or were listed 
bv lHill et all (2005) as having mass values that are unchar- 
acteristic of massive star formation regions. To investigate 
whether statistically significant relationships exist between 
the dust clump properties and the associated methanol 
masers the dust clumps were split into groups depending 
upon the presence or absence of specific associated phenom- 
ena. Four different categorizations were tested, specifically: 

• dust clumps with and without associated 6.7-GHz 
methanol maser emission; 

• dust clumps with 6.7-GHz methanol maser emission 
with and without associated 12.2-GHz methanol maser emis- 
sion; 

• dust clumps with and without detectable radio contin- 
uum; 

• dust clumps associated with more luminous 6.7-GHz 
methanol masers together with those associated with less 
luminous methanol maser sources (sources in this group were 
split into presence/absence data at the median value of 6.7- 
GHz peak luminosity, i.e. presence meaning the 1.2-mm dust 
clump is associated with a 6.7-GHz methanol maser with 
a peak luminosity greater than the median value and vice 
versa for the absence data). 

A Binomial generalize d linear model (GLM) 
(|McCullagh fc Nelderl I1989T ) was then fitted to the maser 
presence/absence data using the dust clump properties as 
predictors for each of the four groups. The significance of 
each individual clump property was tested by fitting all 
possible single term models, and compared by analysis of 
deviance to the null model consisting of only an intercept. 
Stepwise model s election based on the Akaike Information 
Criteria (AIC) l|Burnham fc Anderson! I2002T ) was used 
to select the simplest model with the greatest predictive 
power and is a trade-off between goodness of fit and model 

complexity. 

For each of the 1.2-mm dust clumps (|fflll et al.1 120051 ). 
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G192. 600-0. 048 
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Figure 1. Spectra of all the 12.2-GHz methanol masers detected in the search towards 6.7-GHz methanol masers. 
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Figure 1. - continued 
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Table 2. Absorption features that we detect at 12.2-GHz towards 6.7-GHz methanol masers. The first three columns in this table are as 
in Table [T] Column 4 indicates whether or not emission at 12.2-GHz was detected, followed in columns 5 to 7 by the flux density of the 
absorption, peak velocity (where the peak corresponds to the minimum in the absorption trough) and velocity range of the 12.2-GHz 
absorption features that we detect. 



IVlethcinol mEiser 


RA 


Dec 


12-GHz 


Ads. Flux 


Peak 


Velocity 


(I b) 


(J2000) 


(J2000) 


emission 


Density 


Velocity 


Range 


(degrees) 


(h in s) 


C ' ") 

\ J 


? 


(Jv) 


(km s -1 ) 


(km s~"^) 


G 206. 543-16. 355 2 


05 41 44.15 


-01 54 44.9 


n 


-0.3 


9 


8,11 


G 269. 153-1. 128 1 


09 03 33.46 


-48 28 02.6 


n 


-0.3 


10 


8,12 


G301.136-0.226 1 


12 35 35.14 


-63 02 32.6 


n 


-0.2 


-39 


-42,-37 


G305.200+0.019 1 


13 11 16.93 


-62 45 55.1 


y 


-0.2 


-34 


-40,-26 


Gsos^os+o^oe 1 -' 7 ) 


13 11 13.71 


-62 34 41.4 


y 


-0.4 


-40 


-49,-30 


G318.948-0.196 1 '< 7 ) 


15 00 55.39 


-58 58 52.8 


y 


-0.5 


-37 


-39,-27 


G323.740-0.263 1 << 7 ) 


15 31 45.45 


-56 30 50.1 


y 


-0.6 


-47 


-58,-43 


G331.278-0.188 1 -< 7 ) 


16 11 26.59 


-51 41 56.7 


y 


-0.4 


-77 


-91,-75 


G 000.496+0. 188 1 


17 46 03.96 


-28 24 52.8 


y 


-0.3 


3 


-14,6 


G008.669-0.356 1 


18 06 18.99 


-21 37 32.2 


y 


-0.4 


37 


34,50 


Goos.ess-o.ses 1 ^ 7 ) 


18 06 23.49 


-21 37 10.2 


y 


-0.4 


37 


36,48 


G010.444-0.018 1 '' 7 > 


18 08 44.88 


-19 54 38.3 


y 


-0.3 


70 


63,77 


G010.627-0.384 1 '< 7 ) 


18 10 29.22 


-19 55 41.1 


y 


-0.5 


1 


-5,3 


Gon.gse-o.eie 1 -' 7 ) 


18 14 00.89 


-18 53 26.6 


y 


-0.5 


39 


27,46 


G 012. 025-0. 025 1 ' (7) 


18 12 01.86 


-18 31 55.7 


y 


-0.3 


113 


103,116 


G012.909-0.260 1 '( 7 ) 


18 14 39.53 


-17 52 00.0 


y 


-0.3 


37 


32,47 


G 016. 86-2. 15 3 


18 29 24.411 


-15 16 04.12 


n 


-0.2 


19 


17,21 


G023.437-0.184 1 '< 7 ) 


18 34 39.25 


-08 31 38.5 


y 


-0.7 


101 


90,109 


G023.440-0.182 1 '< 7 ) 


18 34 39.18 


-08 31 25.3 


y 


-0.6 


101 


91,112 


G025.710+0.044 3 '( 8 ) 


18 38 03.148 


-06 24 14.80 


y 


-1.7 


94 


84,100 


G025.826-0.178 3 


18 39 03.631 


-06 24 09.64 


y 


-0.3 


91 


86,95 


G028.305-0.387 3 


18 44 21.988 


-04 17 38.42 


y 


-0.3 


85 


76,87 


G029.956-0.016 3 >( 7 ) 


18 46 03.763 


-02 39 19.69 


y 


-0.4 


100 


90,108 


G029.979-0.047 3 '( 7 ' 


18 46 12.974 


-02 38 58.05 


y 


-0.3 


100 


93,108 


G030.704-0.068 3 


18 47 36.801 


-02 00 48.99 


ii 


-0.4 


91 


89,95 


G030.818-0.057 3 '( 7 ' 


18 47 46.984 


-01 54 19.61 


y 


-0.4 


97 


94,99 


G030.898+0.162 3 


18 47 09.133 


-01 44 10.29 


y 


-0.4 


105 


95,106 


G031.282+0.062 3 -( 7 ) 


18 48 12.385 


-01 26 22.60 


y 


-0.9 


109 


100,116 


G031.411+0.307 3 '< 7 ) 


18 47 34.314 


-01 12 47.14 


y 


-0.2 


99 


96,101 



all of the derived clump properties in each of the four 
groups were tested: the integrated flux density (Jy), peak 
flux density (Jy beam -1 ), source full width at half maximum 
(FWHM) (arcsec), distance (kpc), mass (Mg), radius (pc) 
and H2 number density (cm -3 ). In the following sections we 
consider p- values of less than 0.05 to be statistically signifi- 
cant (i.e. the hypothesis that the single term model provides 
no better fit than the null model consisting only of an in- 
tercept is rejected when the p- value is less than 0.05). Indi- 
vidual p- values, detailed tables and box plots relating to all 
subsequent statistical analysis are presented in Appendix [Al 
A summary of the main results from the statistical analysis 
is presented in Section [3731 



iHill et all l|2009l ) used spectral energy distribution 
(SED) modelling to determine a statistically probable range 
of values f or tempera t ure, m ass and luminosity for around 
half of the IHill et ail ||2005|) sample. We have not incorpo- 



mpli 

rated the data fr omlHill et al. (I2009T) as it is only a subset 



of the sample of IHill et aL n|2005l ) and appears to be bi- 
ased to wards including a large fraction of the more evolved 
sources. IHill et al.l (2009) commented that a large fraction of 
the 1.2-mm dust clumps with no associated methanol masers 
or radio continuum do not have MSX counterparts. As a re- 
sult, a number of these sources were initially included in the 



sample but were later exclud ed as the models were poorly 
constrained (IHill et al.ll2009l). F urthe rmore, 50 of the 1.2- 
mm dust clumps in IHill et alj l|2009l) are associated with 
8- GHz radio continuum emission compared with 71 sources 
in the complete sample of lHill et al] l|2005h . lHill et al.l (|2009h 
find that the masses of sources calculated using the i r SED 
method are similar to those presented in IHill et all (|2005l ) 
with both publications showing that the "MM-only" sources 
are the least massive, followed by the sources associated with 
6.7-GHz methanol masers, and then the sources with radio 
continuum and finally the sources associated with both ra- 
dio continuum and 6.7-GHz methano l maser emission . This 
means that even though the masses in IHill et all (2005) wore 
calculated using a constant temperature of 20 K the calcu- 
lated masses are not si g nifica ntly different from the masses 
presented in IHill et aL I (120091 ) which were calculated using 
temperatures that were determined from their SED mod- 
elling. 

For all dust clump properties with near/far distance am- 
biguities we have used the values calculated at the near dis- 
tance. While it is often considered that, where distance am- 
biguities exist, the near distance is st atistically more likel y 
to be correct th an the far distance (e.g. lSobolev et ai1l2005l ). 
IXu et all (|2009l ) suggest that using the far distances results 
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Figure 1. - continued 



in a better correlation between the locations of methanol 
masers with respect to the spiral arms of the Galaxy. In or- 
der to insure that the assumption that sources are more 
likely to be located at the near distance did not intro- 
duce any significant bias, the distribution of d ust clump 
properties mass and radius of iHill et all ([2005) was com- 
pa red to the distrib u tion o f dust clump prop e rties o bserved 
bv lMookeriea et all (|2004 ). IMookeriea et all |2004) carried 
out 1.2-mm dust clump observations towards the G333 gi- 
ant molecular cloud (GMC) associated wit h RCW106. Thi s 
GMC is located at a distance of 3.6 kpc l|Lockmanlli979l ). 
and all clump properties were calculated at t his distance. 
13 CO observations towards the G333 GMC dBains et all 
2006) showed that, while there was emission detected along 
the line-of-sight, the brightest broadest feature seen in the 
mean velocity profile could be attributed to the target GMC. 



It is therefor e likely that the majori ty of the dust clumps 
observed by IMookeriea et al.l (12004 ) have been correctly 
assigned the distance of the GMC, introducing little bias 
in the clump property distributions. Comparison between 
the distributions of dust cl u mp m ass and radius of the 
d ust clumps from Hill et alj l|2005l ) with the dust clumps 
of IMookeriea etaL ( 2004T ) has shown that the di stribution 
of the clump properties at the near distances of IHill et al.l 
(2005) are similar to the d istribution of the clump properties 
of IMookeriea et all (|2004 ). We therefore expect that our as- 
sumption that the near distances are correct will introduce 
no gross bias into our statistical analysis. Peak luminosities 
of the 6.7-GHz methanol masers were calculated using the 
6.7-GHz methanol maser peak flux density (as integrated 
flux densities are not readily available) and the near kine- 
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Figure 1. - continued 

matic distance of the maser sources as used in Ifiill et ail 
|2005l ). 

In Sections [3T2T21 and EHHI onlv 101 of the 1.2-mm dust 
clumps sources with associated 6.7-GHz methanol masers 
are considered. Five 1.2-mm dust clumps were removed from 
the observed sample of 106 dust clumps with associated 6.7- 
GHz methanol masers for the following reasons; G 259.94- 
0.04 and G 331.28-0.19 have some clump properties miss- 
ing, the 6.7-GHz methanol maser associated with G 284.35- 
0.42 had previously h ad a 12.2-GHz counterpart observed by 
ICaswell et ail (|l995bh that we failed to detect, G 291.58-0.53 
and G 290.374+1.661 were excluded because we detected 
marginal 12.2-GHz emission towards these sources in 2008 
June; but not 2008 December making the se sources i mpos - 
sible to classify and G 6.60-0.08 because Ifflll et all (|2005h 
suggest that the derived mass of this source is uncharacter- 
istic of massive star formation regions and were therefore 
doubtful that it is located at the near distance. 



3.2.1 1.2-mm dust clumps with and without 6.7-GHz 
methanol masers 

Although the observations of iHill et all ((2005) were pri- 
marily targeted towards sites of known 6.7-GHz methanol 
masers, many additional 1.2-mm dust clumps were detected 
in the fields of the observations and have no associated 6.7- 
GHz methanol maser, allowing us to determine if there is any 
statistically significant difference between those dust clump 
sources with and those without associated methanol masers. 
The combined sample size for this analysis is 378. Fits of the 



single term addition Binomial model showed an increasing 
probability of the presence of 6.7-GHz methanol masers as- 
sociated with increasing values of dust clump integrated flux 
density, peak flux density, FWHM, mass and radius. This 
means that any one of them can give an indication of the 
likelihood of having an associated 6.7-GHz methanol maser. 
Results of the single term additions are given in Table [ATI 
The simplest model with the greatest predictive power 
contains the four dust clump properties: integrated flux den- 
sity, peak flux density, FWHM and density. The estimated 
regression relation is 

log—^ — = -2.677 - 0.295xi Iltcgratod + 1.172z Peak 
1 — Pi 

+0.030s; P whm + 0.00000143x Dcnsity 

where reintegrated is the 1.2-mm dust clump integrated flux 
density, Kp ea k is the peak flux density, xfwum is the source 
FWHM, ^Density is the H2 number density and pi is the 
probability of finding a 6.7-GHz methanol maser towards 
the i 1.2-mm dust clump. The regression summary of this 
model is shown in Table lA2l 

However, the misclassification rate for this model is rel- 
atively high. Setting the threshold probability to 0.5 the 
model misclassifies almost two thirds of the 1.2-mm dust 
clumps with associated 6.7-GHz methanol masers as hav- 
ing no associated methanol masers and, while lowering the 
probability threshold allows more of the methanol maser as- 
sociated 1.2-mm dust clumps to be correctly classified, it 
results in a marked increase in the number of 1.2-mm dust 
clumps with no associated methanol masers wrongly classi- 
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fied as having an associated methanol source. A likely ex- 
planation for this is that all of the properties involved in the 
regression model show a significant overlap in the range of 
values for the 1.2-mm dust clump with and without associ- 
ated methanol maser sources (see Fig. IA1|I thereby leading 
to the misclassification of many sources. However , it is also 
possib le that some of the assumptions made by iHill et all 
(2005) in calculating the dust clump properties (like uni- 
form temperature) could account for this. Therefore, a reli- 
able predictive model of 6.7-GHz methanol maser presence 
towards 1.2-mm dust clumps requires more accurately de- 
termined dust clump properties and additional information, 
possibly from other wavelengths. This is in contrast to wa- 
ter maser associat ions with 1.2-mm dust clumps, for which 
iBreen et all |2007) found a promising model that was based 
on only 1.2-mm dust clump radius. 

iBains et al.l (|2009T > co mpared t h e loca tions of 6.7- 
GHz methanol masers l|Ellingsenl 120051 ) with both 
"CO (IBains et all 120061 ) and 1.2-mm dust clumps 
l|Mookeriea et al 2004 ) and used an identical statisti- 
cal analysis to investigate the differences between the 
clum ps with associate d methanol masers and those with- 
out. |Bains^tjjl| (|2009l ) found that 1.2-mm dust clumps with 
associated methanol masers had larger radii, mass and H2 
number densities than those clumps wit h no associat e d 6.7 - 
GHz methanol masers. The results of IBains et all {2009) 
seem to be in contradiction to the results that we find here: 
that there is no statistical difference between the 1.2-mm 
dust clumps H2 number densities for sources with associ- 
ated 6.7-GHz methanol masers compared to those sources 
with no associated 6.7-GHz methanol masers. However, this 
apparent discrepancy can be explained by the considera- 
tion of a few differences between the data sets. C omparison 
of the distr ibution of 1.2-mm dust clumps from IHill et all 
|2005l ) an d iMookeriea et all J2004I ) uncovers the first such 
difference. IMookeriea et all (|2004T ) observed 95 1.2-mm dust 
clumps within a giant molecular cloud (GMC) having H2 
number densities in the range 0.5 to 6 x 10 4 cm -3 . These 
observations had a residual noise limits of 23 mjy for the 
majority of the observed region and up to 50 m jy for the 
northe rnmost section of the GMC. In comparison. IHill et all 
(2005) observed a number of individual sources at a range 
of Galactic longitudes with slightly higher residual noise of 
50 mjy for the majority of sources and up to ~150 mjy 
for some sources. The H2 n umber den s ities d erived for the 
404 1.2-mm dust clumps bv lHill et ail (|2005l ) are all in the 
range 1.4 x 10 3 to 1.9 x 10 6 cm" 3 , a much large range than 
observed by IMookeriea et all (|2004l ). 

IBreen et all (|2007T ) produced a three color GLIMPSE 
image of the GMC and showed that the central section of the 
GMC is dominated by bright mid-infrared emission, trac- 
ing a relatively later evolutionary stage of massive star for- 
ma tion. The majori t y of t he 1.2-mm dust clumps observed 
by IMookeriea et al. I (|2004l ) are also aligned along this sec- 
tion of the GMC, meaning that the maj ority of the 1 . 2-mm 
dust sources will also be more evolved. IBreen et all (|2007i ) 
showed that the methanol masers within the region were 
located towards the periphery of the GMC away from the 
bright mid-infrared emission associated with the central axis 
of star formation, supporting the idea that methanol masers 
tra ce an early evolut ionary stage of massive star formation. 
As IHill et all l|2005l ) targeted many of their observations 



towards methanol masers we expect that the majority of 
sources observed will be at an earlier ev olutionary phase 
than the majority of sources observed by IMookeriea et all 
l(20ul . 

The density of the 1.2-mm dust sources fou nd to be as- 
socia ted with methanol masers in the GMC (|Bains et all 
2009) have densities that would fall in the lower end of 
the range of the H2 number densities we find associated 
with 6.7-GHz m e thano l masers. The small sample size o f 
IMookeriea et all l|2004l) (and therefore IBains et all (|2009l )) 
limits our ability to meaningfully compare the results of 
IBains et all {2009) with what we find here. Furthermore, if 
the difference in the nature of the 1.2-mm ta rget sources of 
IHill et all (|2005l ) and IMookeriea et all l|2004) as well as the 
differing sensitivities of these two sets of observations are 
considered it is not surprising that the results differ; and 
therefore this difference cannot be considered a contradic- 
tion. 



3.2.2 1.2-mm dust clumps and 6.7-GHz methanol maser 
luminosity 

Dust clumps associated with more luminous methanol 
masers and dust clumps associated with 6.7-GHz methanol 
masers with lower luminosities were compared to deter- 
mine if there was any statistically significant dependance 
on 6.7-GHz methanol maser luminosity. The 101 1.2-mm 
dust clumps with appropriate data and associated 6.7-GHz 
methanol masers were split into two groups at the median 
value of 6.7-GHz methanol maser peak luminosity and the 
two groups were compared . In addition to the dust clump 
properties l|Hill et alll2005l ) the significance of the presence 
of an associated OH maser was tested in regions common 
with th e com plete search for OH masers carried out by 
ICaswelll dl998) (this incorporates ~70 per cent of the com- 
plete sample). Fits of the single term addition Binomial 
model to the 1.2-mm dust clump properties showed an in- 
creasing probability of more luminous 6.7-GHz methanol 
masers being associated with dust clump sources that have 
bigger radii, higher mass, lower densities and are located at 
further distances. A summary of the single term addition Bi- 
nomial model can be found in Table lA"3l and is represented 
graphically in the form of box plots in Fig. IA2I 

In order to determine if OH masers are preferentially 
associated with more luminous 6.7-GHz methanol masers 
we then fitted a Binomial model to the presence or absence 
of an OH maser. For this analysis only the sources t h at fel l 
within the range of the complete search of ICaswell (|l998l ) 
were considered. We found that there is a greatly increased 
chance of finding an OH maser associated with the more 
luminous 6.7-GHz methanol masers, the p-value of 0.012, 
corresponding to a confidence level of 98.8 per cent. 

The tendency for the more luminous 6.7-GHz methanol 
masers to be associated with 1.2-mm dust clumps with lower 
H2 number densities was further investigated to see if there 
was an overall trend present in the data. Figure [2] shows the 
log of the 6.7-GHz peak luminosity versus the log of the H2 
number density of the associated 1.2-mm dust clump. We 
have used the technique of linear regression and have found 
that there is a statistically significant negative slope in the 
distribution, with a weak correlation (correlation coefficient 
of 0.46) between the data points. The null hypothesis that 
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Figure 2. Log of 6.7-GHz peak luminosity versus the H2 number 
density of the associated dust clump in cm -3 . The line of best fit 
is also plotted. 

there is zero slope in the data can be rejected (p-value of 
1.16e-06). The line of best fit has the equation 

log(6.7-GHz him) = -0.85[0.16](log(H 2 density)) + 6.25[0.75] 

where '6.7-GHz firm' is the peak luminosity of the 6.7- 
GHz methanol maser and 'H^density' is the H2 number den- 
sity of the associated 1.2-mm dust clumps and the numbers 
contained in the square brackets are the standard errors of 
the slope and the y intercept respectively. It is likely that at 
least some of the scatter present in the distribution of data 
points can be explained by the use of the near kinematic 
distances which are prone to errors. A further contributer 
is li kely to be the variable nature of the mas er sources (see 
e.g. iGoedhart. Gavlard fc van der WaltJl2004D . 

The most parsimonious model produced by fitting a 
GLM to the high and low 6.7-GHz methanol maser lumi- 
nosity associated dust clumps includes both the dust clump 
FWHM and the distance. The estimated regression relation 
is 



= -3.184 + 0.017:tfwhm + 0.489s Dis ta 



Pi 

where ifwhm is the FWHM of the 1.2-mm dust clump and 
distance is the distance to the source. The regression sum- 
mary of this model, including an estimate of errors is shown 
in Table IA41 

The accuracy of this model was tested on the data and 
when the threshold of probability was set to 0.5 the misclas- 
sification rates were found to be promising. Of the 51 1.2- 
mm dust clumps with 6.7-GHz methanol masers with high 
luminosities, the model correctly predicts that 38 of these 
1.2-mm dust clumps will have associated methanol masers 
with high luminosities, but falsely predicts that the remain- 
ing 13 clumps will have associated low luminosity methanol 
masers. Of the 50 1.2-mm dust clumps with associated low 



luminosity methanol masers the model predicts 39 correctly 
and falsely predicts that 11 will have highly luminous as- 
sociated 6.7-GHz methanol masers. It is likely that these 
misclassification rates can be explained by a combination of 
the fact that we are using 6.7-GHz peak luminosity, rather 
than the more appropriate value of 6.7-GHz integrated lumi- 
nosity, the fact that the calculated luminosities are sensitive 
to errors in the adopted kinematic distances and the fact 
that there are a number of sources with similar 6.7-GHz 
peak luminosities around the median value. 



3.2.3 1.2-mm dust clumps with associated 6.7-GHz 
methanol masers with and without 12.2-GHz 
methanol masers 

The properties of the 1.2-mm dust clumps associated only 
with 6.7-GHz methanol masers were compared to 1.2- 
mm dust clumps associated with both 6.7- and 12.2-GHz 
methanol masers; the sample size for this analysis was 101 
dust clumps (those that were excluded are shown in col- 
umn 12 of Table [1] and discussed in Section I3.2.2|l . Fits of 
the single term Binomial model show that there is a sta- 
tistically significant increase in the likelihood of finding a 
12.2-GHz methanol maser associated with increasing 6.7- 
GHz methanol maser flux density, 6.7-GHz methanol maser 
peak luminosity and decreasing values of 1.2-mm dust clump 
density. The results of the single term additions are in Ta- 
ble IA5I The significance of the presence of an OH maser 
associated with the 6.7-GHz methanol masers was tested 
on the s ample within the search region for OH masers of 
ICaswelll (|1998i ) (this incorporates ~70 per cent of the com- 
plete sample). This resulted in a p-value of 0.025, meaning 
that 12.2-GHz methanol masers are preferentially found to- 
wards 6.7-GHz methanol masers with associated OH masers. 

The luminosity of the 12.2-GHz methanol masers were 
investigated to determine if individual source luminosity had 
any relationship with either the association with 8-GHz ra- 
dio continuum or the luminosity of the associated radio con- 
tinuum. We found that there is no relationship between ei- 
ther the presence of an associated radio continuum source 
or the luminosity of associated radio continuum and the lu- 
minosity of the associated 12.2-GHz methanol maser. 

The luminosities of the 12.2-GHz methanol masers were 
compared to the H2 number densities of the 1.2-mm dust 
clump using the same method as for 6.7-GHz peak lumi- 
nosity in Section 13.2.21 Figure [3] shows a plot of the data 
overlaid with the line of best fit. We find that there is a very 
weak correlation in this data set, with a correlation coeffi- 
cient of 0.25. Similarly to Fig. [2] we expect that the largest 
contributions to the scatter in Fig. [3] are the use of the kine- 
matic distances and maser variability. However, there is an 
overall negative slope in the data that is statistically signif- 
icant, the p-value is 0.05 meaning that the null hypothesis 
that there is zero slope can be rejected. The line of best fit 
has the equation 



log(12.2-GHz lum) = -0.56[0.28](log(H 2 density)) + 4.16[1.23], 

where '12.2-GHz lum' is the integrated luminosity of 
the 12.2-GHz methanol maser and the '^density' is the 
H2 number density of the associated 1.2-mm dust clumps 
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Figure 3. Log of 12.2-GHz integrated luminosity versus the H2 
number density of the associated dust clump in cm -3 . The line 
of best fit is also plotted. 



and the numbers contained in the square brackets are the 
standard errors of the slope and the y intercept respectively. 

The most parsimonious model produced by fitting a 
GLM to the 12. 2- GHz methanol maser presence/absence 
data, includes only the value of the 6.7-GHz peak flux den- 
sity. The estimated regression relation is 



lo 



>grr^— = -0.41058 + 0.03571x 6 . 7 ghzfiu X , 
I- Pi 



where 3%.7GHzFiux is the flux density of the 6.7-GHz methanol 
maser peak in Jy. The regression summary of this model, 
including an estimate of errors is shown in Table IA6I The 
physical implication of this model is that 6.7-GHz methanol 
masers with peak flux densities of ~11.5 Jy have a probabil- 
ity of 0.5 of having a 12. 2- GHz methanol maser counterpart 
with a flux density >0.55 Jy, with the probability increas- 
ing with increasing values of 6.7-GHz methanol maser flux 
density and decreasing for lower values. 

For the 63 dust sources that do have associated 6.7- and 
12.2-GHz methanol masers the model predicts (when the 
probability threshold is set to 0.5) that 44 of the 6.7-GHz 
methanol masers will have associated 12.2-GHz methanol 
maser emission and wrongly predicts that the remaining 19 
will not have associated 12.2-GHz methanol maser emission. 
Applying the model to the 38 6.7-GHz methanol masers with 
no detectable 12.2-GHz emission, the model predicts 29 of 
these correctly and falsely expects the remaining 9 to have 
associated 12.2-GHz methanol maser emission. 

This analysis was then repeated after limiting the sam- 
ple to sources where the associated 6.7-GHz methanol maser 
was stronger than 9 Jy (more than 16 times our detection 
limit of 0.55 Jy). This was done to determine if the results of 
the analysis of the full sample was due to a sensitivity bias 
that may have been introduced by the fact that 12.2-GHz 
methanol masers are usually weaker th an their associated 
6.7-GHz emission. ICaswell et aT] (|l995bh found that the me- 



dian ratio of 12.2- to 6.7-GHz methanol maser intensity was 
1:9 (taking into account 12.2-GHz non-detections) which we 
have doubled and adopted as the 6.7-GHz methanol maser 
cut off threshold. Therefore, methanol masers with peak flux 
densities of less than 9 Jy were removed from the sample, 
leaving 61 dust clumps, and the statistics were recalculated. 
The removal of the weaker 6.7-GHz methanol maser sources 
had little effect on the results of the statistics with the 1.2- 
mm dust clump density, 6.7-GHz methanol maser peak lumi- 
nosity and flux density remaining the only properties show- 
ing statistically significant differences. 



3.2.4 Statistical analysis of 1.2-mm dust clumps with and 
without 8- GHz radio continuum 

In order to determine if there was any difference between 
the properties of 1.2-mm dust clumps with and without de- 
tectable radio continuum we also fitted a Binomial GLM 
to the radio continuum presence/absence data. 378 1.2-mm 
dust clumps were used in t his analysis, 71 of which have 
detectable radio continuum l|Hill et al.ll2005l . and references 
therein). Fits of the single term Binomial model showed that 
all tested 1.2-mm dust clump properties, with the excep- 
tion of distance, were statistically significant. In the case 
of 1.2-mm integrated and peak flux densities, radius, mass 
and FWHM, increasing values of the dust clump proper- 
ties showed a greater likelihood of being associated with 
detectable radio continuum, whereas in the case of the H2 
number density it is the smaller values that show an in- 
creased likelihood of associated radio continuum. Results of 
the single term additions are in Table IA8I and that data is 
shown in the form of box plots in Fig. IA5I 



3.3 Summary of statistical results 

Here we briefly summarize the main results from the statis- 
tical analysis carried out in the previous section, Section[IO] 



• 1.2-mm dust clumps with associated methanol masers 
have higher values of integrated and peak flux densities, 
mass and radius than those with no associated 6.7-GHz 
methanol maser. 

• 1.2-mm dust clumps with associated 6.7-GHz methanol 
masers with high peak luminosities are more likely to be 
associated with OH masers and have larger values of mass 
and radius and lower values of density than those clumps 
that are associated with 6.7-GHz methanol masers with low 
luminosities. 

• 1.2-mm dust clumps associated with both 6.7- and 
12.2-GHz methanol masers are more likely to be associated 
with OH masers and have lower densities than 1.2-mm dust 
clumps with just 6.7-GHz methanol maser associations. 

• The luminosity of 12.2-GHz methanol masers is not de- 
pendent on either the presence of associated radio contin- 
uum, or the luminosity of associated radio continuum. 

• 6.7-GHz methanol masers with 12.2-GHz methanol 
maser counterparts have larger flux densities and higher 
peak luminosities than 6.7-GHz methanol masers with no 
12.2-GHz counterpart and this is not likely to be a result of 
a sensitivity bias. 

• 1.2-mm dust clumps with associated radio continuum 
have higher values of mass, radius, peak and integrated flux 
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density and lower values of density than clumps without 
detectable radio continuum. 



4 DISCUSSION 

4.1 Comparison between 6.7- and 12.2-GHz 
methanol masers 

We have made a preliminary investigation of the relative 
intensities and velocity ranges of the 6.7- and 12.2-GHz 
methanol maser sources within this sample. Detailed analy- 
sis of these properties, as well as spectral structure, will be 
carried out in a subsequent paper which will present 12.2- 
GHz observations towards all 6.7-GHz methanol masers de- 
tected in the southern hemisphere comp onent of the Parke s 
Methanol Multibeam (MMB) Survey (|Green et al J [2009V 
The MMB survey for 6.7-GHz methanol masers is the most 
sensitive survey yet undertaken for young high-mass stars in 
the Galaxy and is complete within 2 degrees of the Galactic 
plane. 

Investigation of the peak flux densities of the 6.7- and 
12.2-GHz methanol maser associations reveals that there 
are three instances where the flux density of the 12.2-GHz 
methanol maser is stronger than the reported values of the 
associated 6.7-GHz methanol maser emission. Given that 
interstellar masers are intrinsically variable it is likely that 
in at least some of these instances the 6.7-GHz methanol 
masers are now stronger than their reported values from 
several years prior to the 12.2-GHz methanol maser obser- 
vations. 

We have made a comparison of the 6.7- and 12.2-GHz 
methanol maser peak flux densities (where emission was de- 
tected at both frequencies) and find that the median ratio 
of 12.2- to 6.7-GHz methanol flux density is 1:5.9 with an 
average ratio of 1:3.1. We additionally find that the most 
extreme ratios of 12.2- to 6.7-GHz methanol maser peak 
flux density are 1:180 and 1:0.56, highlighting the fact that 
when the 12.2-GHz emis sion is stronger it is on ly marginally 
stronger. In comparison, ICaswell et all (|l995bh found a me- 
dian ration of 12.2- to 6.7-GHz methanol maser peak flux 
density of 1:5.4. 

We have compared the ratio of 6.7- to 12.2-GHz peak 
flux density to both the peak luminosity of the 6.7-GHz 
methanol maser and the integrated luminosity of the 12.2- 
GHz methanol masers. We find that there is a general trend 
whereby the 6.7-GHz methanol masers with high peak lu- 
minosities have higher ratios of 6.7- to 12.2-GHz methanol 
maser flux density (i.e. the more luminous 6.7-GHz methanol 
masers have relatively weaker emission at 12.2-GHz). Com- 
parison between source flux density ratios and 12.2-GHz 
methanol maser luminosities calculated from the integrated 
flux densities of the 12.2-GHz methanol maser sources show 
the complementary trend (i.e. the less luminous 12.2-GHz 
methanol masers have relatively stronger emission at 6.7- 
GHz). 

Plots of the flux density ratios versus 6.7-GHz peak lu- 
minosity and 12.2-GHz integrated luminosity are shown in 
Figs [4] and [5] to demonstrate our results quantitatively. We 
have used linear regression to show that the trends in the 
data are statistically significant. For the 6.7-GHz peak lu- 
minosity, we find a slope of 0.25 with a standard error of 



0.08 and a p- value of 0.003 which allows us to reject the null 
hypothesis that there is zero slope. The correlation coeffi- 
cient between the points is 0.36 which is classified as a weak 
correlation. For the 12.2-GHz integrated luminosity we find 
a slope of -0.15 with a standard error of 0.07 and a p- value 
of 0.03 which similarly means that the slope of the line is 
statistically significant and non-zero. The correlation coef- 
ficient is 0.27 which also corresponds to a weak correlation 
between the data points. The significant scatter present in 
both data sets are likely to be due to a combination of a 
few factors. The first is that the calculated luminosities are 
heavily dependent on the derived kinematic distances for 
the sources, which in many instances have large ambigui- 
ties. We have assumed the near distances for sources here 
as for all previous analysis which, for some sources, may in- 
troduce large errors into the calculated source luminosities. 
It is also likely that the 6.7-GHz methanol maser flux den- 
sities reported in the literature would not correspond to the 
true flux density of these sources if measured at the time 
of the 12.2-GHz methanol maser observations. We therefore 
expect that source variability also contributes to the scatter 
in Figs [4] and [5] It is expected that these influences will have 
an effect on our data; however, we would expect that they 
would act to weaken trends not acting in a systematic way 
to strengthen them. It is therefore likely that accurate dis- 
tance measurements and single epoch observations of both 
methanol maser transitions would result in a more tightly 
correlated relationship, but, this is something for future in- 
vestigation. 

In section 3.2.2 we showed that the 6.7-GHz peak lu- 
minosity increases with decreasing dust clump density. We 
also show that centimetre radio continuum emission is pref- 
erentially found towards lower density dust clumps, which 
suggests that the dust clump density decreases as the source 
evolves (see section 4.3 for a more detailed discussion of our 
results and their implications for an evolutionary sequence 
for high-mass star formation). These results (and others) 
combined suggest that the 6.7-GHz methanol peak luminos- 
ity increases as the high-mass star formation it is associated 
with evolves. That the average ratio of 6.7- to 12.2-GHz 
peak flux density increases as the 6.7-GHz peak luminos- 
ity increases implies that the rate of increase of the 12.2- 
GHz peak flux density as the source evolves is less than that 
for the 6.7-GHz transition. Modeling of methanol masers by 
Cragg et al. (2005) shows that the 6.7-GHz transition is es- 
sentially always stronger than the 12.2-GHz transition and 
our results suggest that as the underlying star formation re- 
gion evolves the relative intensity of the 6.7-GHz transition 
increases more quickly. In particular, Figures 2 & 3 of Cragg 
et al. (2005) show that the intensity of the 12.2-GHz transi- 
tion decreases more rapidly with increasing gas temperature 
than the 6.7-GHz transition does, which is consistent with 
an evolutionary trend for the 6.7- to 12.2-GHz flux density 
ratio to increase. 

Inspection of the velocity ranges of the 6.7-GHz 
methanol masers as well as their associated 12.2-GHz 
methanol maser emission has shown that of the 68 12.2- 
GHz methanol masers that we observed, 49 sources show 
the peak of their emission at the same velocity as the 6.7- 
GHz methanol maser, to within the spectral resolutions of 
the respective observations. For the remaining 16 sources we 
find the 12.2-GHz peak velocity within the velocity range 
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Figure 4. Log of the 6.7-GHz methanol maser peak flux density 
divided by the 12.2-GHz methanol maser peak flux density versus 
the log of the 6.7-GHz peak luminosity, overlaid with the line of 
best fit. 
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Figure 5. Log of the 6.7-GHz methanol maser peak flux density 
divided by the 12.2-GHz methanol maser peak flux density versus 
the log of the 12.2-GHz integrated luminosity, overlaid with the 
line of best fit. 



of the detectable 6.7-GHz methanol maser emission and is 
coincident with a secondary 6.7-GHz methanol maser fea- 
ture. In the majority of cases we find that the 12.2-GHz 
methanol masers have narrower velocity ranges than the 
6.7-GHz methanol maser counterparts. The average velocity 
range of the 6.7-GHz methanol masers with associated 12.2- 
GHz methanol masers is 9 km s" 1 and the average velocity 
range of the 12.2-GHz methanol masers that we observe is 
3.4 km s _1 . The median value of 12.2-GHz velocity range of 



the sources that we observed is 5.4 km s 1 which is remark- 
ab ly similar to the media n value of the 131 sources reported 
bv lCaswell et~ai1 (|l995bh as 5 km s 



4.2 Absorption at 12.2-GHz 

29 of our target 6.7-GHz methanol masers show absorp- 
tion at 12.2-GHz, of w hich 24 also show 12.2 -GHz maser 
emission. In comparison ICaswell et all (|l995bh detected 19 
absorption spectra toward 238 6.7-GHz methanol masers. 
We suspect that absorption accompanying emission in some 
cases may be an artifact of out processing method (see Sec- 
tion [2]). 



4.3 Evolutionary sequence for sources 

The relatively large number of class II methanol maser 
sources in our Galaxy (esimated to be of order 1000) , means 
that they must be associated with objects covering a range 
of stellar masses and also have an appreciable "lifetime" . It is 
generally assumed that the properties of the masers and the 
complementary data at different wavelengths will depend 
upon both the age of the source and its mass. At present 
we do not have sufficient information to disentangle these 
two factors. Circumstantial evidence can sometimes be used 
to infer one or the other as the primary factor. However, 
some of the attributes ascribed to evolutionary factors may 
be partially, or ev e n larg ely mass related (and vice versa). 

iPurcell et alj l|2009l ) propose that the difference be- 
tween radio-loud and radio-quiet massive young stellar ob- 
jects (YSOs) may be attributed to evolution, with the radio- 
loud sources corresponding to a later stage of evolution. Our 
statistical comparison between 1.2-mm dust clumps with as- 
sociated radio continuum from ultra-compact (UC) Hn re- 
gions at 8 GHz and those without has shown that 1.2-mm 
dust clumps with detectable radio continuum in general have 
bigger radii, masses, flux densities as well as lower H2 num- 
ber density. We therefore propose that these 1.2-mm dust 
clump properties give an indication of comparative stages of 
evolution (i.e. sources with larger radii, masses, flux densi- 
ties and smaller H2 number densities are more likely to be 
at a later stage of evolution than sources which have smaller 
value s of radii, mas s , flux density and larger values of den- 
sity). IPurcell et all (|2009l ) similarly found that YSOs with 
associated 8-GHz radio continuum are more massive than 
those without detectable radio continuum. 

The evolutionary interpretation of the results of the 
statistics is contingent on the assumption that a large num- 
ber of the 1.2-mm dust clump sources with no detectable 
radio continuum will go on to form massive stars. While 
it is likely that some of the 1.2-mm dust clumps with- 
out associated radio continuum are 'failed' cores and will 
not form massive stars, we know that a large number of 
them will do so, since they exhibit other tracers that are in- 
dicative of ma s sive star formation such as methanol masers 
(ICaswell|[2009l; iMinier et al.ll2003l: TWalsh et al.l Il998l . Il997l ; 
ICaswell et alj Il995al; Szvmczak et al.l I2002T ) and ammonia 
( Longmore et al.ll2007t ). lHill et al.1 d2009h concl uded from in- 
vestigation of a sub-set of the Irlill et alj l|2005h sample that 
the 1.2-mm dust clumps with no overt signs of massive star 
formation displayed similar characteristics to those 1.2-mm 



Quantifying high-mass star formation evolutionary schemes 23 



dust clumps associated with known massive star formation 
regions. Therefore, as the number of sources in this sample 
is large we believe that the possibility of the inclusion of 
some 'failed' cores is inevitable, but the number will be few 
and so will have little effect on our statistics. 

lLongmore et al.l (|2007T ) detected 24-GHz continuum 
emission towards some of the 1.2-mm dust clump sources 
where no 8 - GHz continuum emission was detected by 
IWalsh et all l|l99St ). It is therefore likely that some of the 
sources that we have classified as having no associated con- 
tinuum emission are in fact associated with radio continuum 
from hyper-compact (HC) Hn regions which are often op- 
tically thick at centimetre wavelengths and are associated 
with an earlier evolutionary phase of massive star fo r matio n 
than UCHu regions. According to lLongmore et all l|2007h . 
this means that sources with detectable radio continuum at 
24-GHz and not 8-GHz are part of another class of object 
that is more evolved than sources with no detectable radio 
continuum and less evolved than sources with detectable 
continuum at centimetre wavelengths. This distinction also 
means that we have not introduced any bias into our anal- 
ysis by not accounting for sources that are in the HCHn 
region stage. 

Interpretation of the relative evolution of sources as- 
sociated with the different maser properties can be made 
by comparing the 1.2-mm dust clump properties associated 
with these sources to the results of the analysis with UCHn 
regions. We find that 1.2-mm dust clumps with associated 
6.7-GHz methanol masers have higher flux densities, radii 
and mass than those with no associated 6.7-GHz methanol 
maser. This suggests that 6.7-GHz methanol masers are 
present prior to any significant change in source density that 
occurs along side the formation of UCHn regions. H o wever , 
the 1.2-mm dust clump observations of iHill et all (|2005l ) 
have a relatively low resolution (24") and so only reflect the 
large scale structures. Hence significant changes in source 
structure and physical conditions are likely to be occurring 
on smaller spatial and temporal scales, but these will only 
be revealed by higher resolution observations of the dust 
properties. 

Our analysis of 6.7-GHz methanol maser peak luminos- 
ity and 6.7-GHz methanol masers associated with 1.2-mm 
dust clumps shows that methanol masers with high peak 
luminosities are more likely to be associated with 1.2-mm 
dust clumps of larger mass and clump radii and lower H2 
number densities than 6.7-GHz methanol masers with lower 
peak luminosities. While we find that there is no statistical 
difference between the flux densities of 1.2-mm dust clump 
sources with varying levels of 6.7-GHz methanol maser peak 
luminosity, we find that the properties of the 1.2-mm dust 
clumps associated with 6.7-GHz methanol masers with high 
luminosities are similar to those 1.2-mm dust clumps that 
are associated with a later stage in evolution. We addi- 
tionally find that more luminous 6.7-GHz methanol masers 
are more likely to be associated with OH masers which 
are known to trace a generally later evolutionary phase of 
massive star formation than isolated m ethanol masers (e.g. 
iForster fc Caswelll Il989l : ICaswelll 1 19971 ). Therefore we ex- 
pect that there is a trend between 6.7-GHz methanol maser 
peak luminosity and evolutionary phase, with the more lumi- 
nous 6.7-GHz methanol maser peaks tending towards a later 
phase in evolution than 6.7-GHz methanol masers with low 



peak luminosities. This important result will deserve further 
investigation using luminosities that have been calculated 
from integrated maser flux densities. 

We find that 1.2-mm dust clump sources associated 
with both 6.7- and 12.2-GHz methanol masers are more 
likely to have lower densities than the sources which 
are associated with only the 6.7-GHz maser transition 
of methanol. Pumping models of methanol masers (e.g. 
ICragg. Sobolev fc Godfrey! l2005h predict that the 6.7- and 
12.2-GHz transitions are inverted over very similar ranges 
of parameter space. The absence of 12.2-GHz masers as- 
sociated with a large fraction of 6.7-GHz sources (40 %) 
implies that the physical conditions are commonly close 
to the point where the 12 .2 -GHz masers switch on or off. 
ICragg. Sobolev fc Godfrey! (|2005l ) also find that 6.7-GHz 
methanol masers can exist at higher densities than 12.2- 
GHz methanol masers, consistent with the results of our 
analysis. We also find that 12.2-GHz methanol masers are 
more likely to be associated with 6.7-GHz methanol masers 
with higher flux densities and peak luminosities. This re- 
sult is not surprising as all previous observations of both 
12.2- and 6.7-GHz methanol masers have shown that the 
12.2-GHz emission is rarely brighter than the associated 
6.7-GHz methanol maser emission. However, when we less- 
ened the chance of a sensitivity bias by removing all of the 
weak 6.7-GHz methanol masers from our analysis, we find 
that our results do not change. This means that there is 
likely to be a delay between when the 6.7-GHz methanol 
maser 'switches on' and when the 12.2-GHz methanol maser 
'switches on'. Furthermore, 12.2-GHz methanol masers are 
more likely to be found toward 6.7-GHz methanol masers 
with associated OH masers. We therefore suggest the 12.2- 
GHz methanol masers are present sometime after the onset 
of 6.7-GHz methanol maser emission. 

Ivan der Waltl ((2005) estimates the lifetime of 6.7-GHz 
methanol masers to be in the range 2.5 x 10 4 to 4.5 x 
10 4 years, combining this with our detection rate of asso- 
ciated 12.2-GHz methanol masers we expect that the cor- 
responding lifetime of 12.2-GHz methanol masers will be 
in the range 1.5 x 10 4 to 2.7 x 10 4 years. This inference 
assumes that 12.2-GHz methanol masers are never with- 
out associated 6.7-GHz methanol maser emission. However, 
given that previous searches have not resulted in serendip- 
itous detections of 12.2-GHz methanol masers without ac- 
companying 6.7-GHz methanol maser emissio n and the fact 
that the current search and previ ous searches l|Caswell et alj 
ll995bl ; lBlaszkiewicz fc Kus!l2004l ') have shown that 12.2-GHz 
methanol maser emission is only rarely brighter that the as- 
sociated 6.7-GHz maser emission, this is not an unreasonable 
assumption. 

While our simplified evolutionary results appear to be 
self consistent and reliable the full picture is likely to be 
more complex. We found that 1.2-mm dust clumps with 
detectable radio continuum have larger radii, masses, flux 
densities and smaller H2 number densities and suggest that 
these properties can be used to give an indication of the 
probable evolutionary stage. However, it is likely that some 
of the more massive dust clumps will go on to form higher 
mass stars and are therefore more likely to have detectable 
radio continuum earlier than some of the less massive cores. 
Furthermore the flux density of the radio continuum will 
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depend on the dust fraction and the type of star that the 
1.2-mm dust clump is associated with. 

The H2 number densities of the 1.2-mm dust clumps are 
central to a number of our evolutionary arguments. There 
are several factors that directly affect the accurate determi- 
nation of this quantity. Both the 1.2-mm dust clump mass 
and radius rely heavily on the accuracy of the distance mea- 
surements. However, it would be unlikely that the suscepti- 
bility of the radius and mass to errors would have a system- 
atic effect that resulted in the strengthening of a relation- 
ship between radius and some other property, such as maser 
presence - it would be much more likely to dilute such a 
relationship. We therefore believe that the only effect that 
potential errors in the determined 1.2-mm dust clump radius 
(and therefore the H2 number density) would have on our 
statistics is to weaken true relationships rather than falsely 
cau se correlations be tween properties. 

iHill et al.l (|2005h adopted a constant tempera te of 20 K 
to calc ulate the masses of the 1.2-mm dust clumps. [Hill et all 
(2009) have shown that in the vast majority of cases this is 
an overestimate but t hat the distribut ion of masses is sim- 
ilar to that found in IHill et al.l ((2005). As we expect less 
evolved 1.2-mm dust clumps to be colder than the more 
evolved sources that show detectable 8-GHz radio contin- 
uum emission, errors introduced by the assumption of a con- 
stant temperature would result in the less evolved sources 
having underestimated masses relative to the more evolved 
sources. Therefore, the use of accurate temperature mea- 
surements would strengthen our finding that more evolved 
sources have lower H2 number densities. Due to this depen- 
dance however, we are unable to separate the effects that 
the H2 number density and temperature have on the evolu- 
tionary phase and expect that temperature would also play 
a large role in indicating the relative evolutionary phase of 
the sou rces. 

As lHill et al.l (|2005l ) targeted their observations towards 
sources suspected of undergoing massive star formation, 
there is a possibility that an unbiased search would result in 
a distribution of source properties that differs from this sam- 
ple. However, as discussed in Section f3. 2. II the distribution 
range of t he 1.2-mm dust cl ump properties from the targeted 
sample of IHill et al.l (|2005h greatly exceed s the distribution 
range from the unbiased search of a GMC l|Mookeriea et all 
120041 '). Furthermore, as the number of serendipitous detec- 
tions of 1.2-mm dust clumps within the fields of targeted 
objects is large (253), the likelihood of biases being intro- 
duced by the targeted nature of this sample is low. 

Another point that needs to be considered is that a dust 
clump is unlikely to form only one star. Some of the confu- 
sion in our predictive models can be explained by the fact 
that often, there are clusters of stars forming within the one 
dust clump, at different evolutionary stages and with differ- 
ent mass ranges. This fact is evident in our sample because 
there are multiple instances where there is more than one 
6.7-GHz methanol maser within the one dust clump. 

Combining all of our statistical analyses we find that as 
the 1.2-mm dust clump becomes more massive (n.b. this is 
mostly likely a consequence of the relatively underestimated 
masses of the less evolved sources, a consequence of the con- 
stant temperature assumption, rather than a real increase in 
mass as discussed above), and increases in radius and flux 
density, a weak 6.7-GHz methanol maser forms. Then, as 



the clump becomes less dense and continues to increase in 
mass, radius and flux density, the 6.7-GHz methanol maser 
increases in flux density and the 12.2-GHz methanol maser 
forms. At this stage the detection of associated OH masers 
and 8-GHz radio continuum becomes more likely. We there- 
fore propose the following scenario for evolutionary sequence 
of sources from less evolved to more evolved: (i) 1.2-mm 
dust clumps with no associated methanol maser emission or 
8-GHz radio continuum, (ii) 1.2-mm dust clumps with asso- 
ciated 6.7-GHz methanol masers, (iii) 1.2-mm dust clumps 
with associated methanol maser emission at both 6.7- and 
12.2-GHz. 8-GHz radio continuum is observed toward 29 of 
our sample of 113 methanol masers and we find that 19 of 
these also have associated 12.2-GHz methanol maser emis- 
sion (this is not statistically significant p-value=0.24 from 
a chi-sqared test). We therefore postulate that 8-GHz radio 
continuum generally becomes detectable toward the end of 
the 6.7-GHz methanol maser only stage, prior to the onset 
of 12.2-GHz methanol maser emission if the mass of the star 
is sufficiently large. 

In Fig. [|J] we present our proposed evolutionary se- 
quence for masers in massive star formation r egions. We have 
used t he evolutionary sequence presented in lEllingsen et all 
(2007) as a starting point, adjusted it according to our new 
data, and placed quantitative estimates on the relative life- 
times of the different speci es using additional maser obser- 
vations (e.g. [Caswell 2003). The estimates of the lifetimes 
of the different maser species have bee n extrapolated from 
the 6.7-GHz methanol maser lifetime of lvan der Waltl (120051 ) 
and detection statistics for the different maser species. 

4-3.1 6.7- & 12.2-GHz methanol maser association with 
radio continuum 

There are 29 1.2-mm dust clumps that are associated with 
both methanol maser emission and 8-GHz radio continuum 
within our sample. 19 of these sources are associated with 
both 6.7- and 12.2-GHz methanol maser emission while the 
remaining 10 are associated with only 6.7-GHz methanol 
masers. We have compared these two groups of sources in 
order to determine if there are any differences in the 1.2- 
mm dust clump properties. We find that the clumps with 
both 12.2 and 6.7-GHz methanol masers have large FWHM 
(p- value 0.01), lower densities (p- value 0.07) and higher 6.7- 
GHz methanol maser peak luminosity and peak flux den- 
sities. We have also investigated the 8-GHz radio contin- 
uu m luminositi e s (cal culated from the flux densities given 
m IWalsh et all l|l998h ) and found that there is no differ- 
ence between either the flux density or the luminosity of the 
continuum between those sources with associated 12.2-GHz 
methanol masers and those without. This implies that the 
presence of an associated 12.2-GHz maser is much more re- 
liant on suitable densities rather than either the presence or 
strength of an associated radio continuum source. In general, 
however, 1.2-mm dust clumps with associated radio contin- 
uum emission are more likely to satisfy the density ranges 
for emission from the 12.2-GHz methanol maser transition. 

4.4 Association with GLIMPSE point sources 

Positions of the 6.7-GHz methanol masers were compared 
to the positions of infrared point sources listed in the 
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Figure 6. Evolutionary sequence for masers associated with massive star formation regions. Lifetimes represented by dotted lines indicate 
that they have been estimated from the literature along with unpublished data. The lifetime of the UCHn region is represented by a 
dashed line as the onset of radio continuum emission is particularly dependent on the mass of the associated star. The double arrow 
head on the far end of the UCHn region lifetime is to show that it will persist past the end of the range of the plot. 



GLIMPSE point source catalogue as well as the GLIMPSE 
archive. As most methanol maser positions are known to 
a high precision we assume an association between the 
GLIMPSE sources and the methanol masers if they lie 
within 2 arcsec of each other. We find that 90 of the 
methanol masers that we observe lie within the regions cov- 
ered by GLIMPSE of which 51 (57 %) are associated with ei- 
ther a GL I MPSE point source or GLIMPSE archive source. 
Ellingscn (2006j) found 68 % of a sample of 56 methanol 
masers had associated GLIMPSE sources. The difference in 
detection rates is not statistically significant. 

Fig. [7] shows a plot of the [3.6]-[4.5] /jm versus [5.8]- 
[8.0] /xm colours of the GLIMPSE point sources associated 
with the methanol masers in our sample. For sources to 
be present on this plot they needed to have flux density 
measurements for all four of the IRAC bands, limiting the 
number of 6.7-GHz methanol maser associated GLIMPSE 
sources in the plot to 13, 8 of which are also associ ated with 
12.2-G Hz methanol maser emission. Similarly to Ellin gsen 
(2006), we find that the GLIMPSE sources associated with 
the methanol masers lie above the majority of the compari- 
son sources in our colour-colour plot. 

We have compared the colours of the GLIMPSE sources 
associated with both 6.7- and 12.2-GHz methanol masers 
with those associated with only 6.7-GHz methanol masers 
and find that there is no statistically significant difference 
between them. We will not be able to determine if this is 
a consequence of our small sample size or rather a real ef- 
fect, until we can compare the colours for a larger number 
of sources. However, if our results are representative then 
this implies that the masers themselves are much more sen- 
sitive to evolutionary changes than is the associated infrared 
emission of wavelengths < 10 /im. 



5 CONCLUSION 

We have searched 113 known 6.7-GHz methanol masers 
for associated 12.2-GHz methanol maser emission with the 
ATNF Parkes 64-m radio telescope. We have detected 12.2- 
GHz methanol masers towards 68 of these sources; a detec- 
tion rate of 60 percent. We estimate the life-time of 12.2- 
GHz methanol maser emission to be of the order of 2.1 x 
10 4 years. We find that for the majority of sources, the peak 
velocity of the 12.2-GHz methanol maser emission is the 
same as the 6.7-GHz methanol maser and, for the remainder, 
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Figure 7. Colour-colour plot of GLIMPSE point source data. 
Methanol maser sources showing both 6.7- and 12.2-GHz emission 
are represented by blue squares and sources with only the 6.7-GHz 
methanol maser transition are represented by the red circles. The 
black dots represent all of the GLIMPSE point sources within 30 
arcmin of I = 326? 5, b = 0?0. 

the 12.2-GHz peak emission is coincident with a secondary 
feature present in the 6.7-GHz methanol maser spectrum. 
Comparison between the peak flux densities of the respec- 
tive methanol maser transitions has shown that median ratio 
of 6.7- to 12.2-GHz peak flux density is 1:5.9, similar to pre- 
vious searches. We additionally find that the ratio of 6.7- 
to 12.2-GHz methanol maser peak flux density has some 
dependence on 6.7-GHz methanol maser luminosity, with 
the more luminous 6.7-GHz sources having comparatively 
weaker 12.2-GHz methanol maser emission and vice versa. 

All of the target sources have prev iously been obse rved 
at 1.2- mm for dust continuum emission ilHill et al.ll2005l ). We 
have carried out statistical analysis of the dust clump prop- 
erties of sources in the following four subgroups compared 
the 1.2-mm dust clump properties: dust clumps with and 
without associated 6.7-GHz methanol maser emission; dust 
clumps associated with highly luminous 6.7-GHz methanol 
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masers compared with those of low luminosity; dust clumps 
associated with 6.7-GHz methanol maser emission with and 
without associated 12.2-GHz methanol maser emission; and 
dust clumps with and without associated radio continuum. 
Our statistics show that there are many differences between 
the sources that fall into these subgroups; consistent with 
the different source combinations and characteristics being 
associated with different stages of massive star formation. 
An evolutionary scenario for the common maser species as- 
sociated with high-mass star formation regions is presented. 

The positions of the 6.7-GHz methanol masers have 
been compared to sources in the GLIMPSE point source cat- 
alogue. We find that 57 per cent of the 6.7-GHz methanol 
masers that lie within the GLIMPSE regions are associ- 
ated with a GLIMPSE point source. Comparison between 
the GLIMPSE colours associated with sources that exhibit 
emission at both 6.7- and 12.2-GHz and those exclusively 
at 6.7-GHz has shown that there is little difference between 
the colours of the associated point sources. 

We will present detailed investigation of the relative in- 
tensities and velocity ranges of 12.2-GHz methanol masers 
detected towards the 6.7-GHz methanol masers detected 
in t he southern hemis phere component of the MMB sur- 
vey dGreen et alj|2009h in a subsequent paper. The near si- 
multaneous observations of hundreds of 6.7- and 12.2-GHz 
methanol masers collected from the MMB survey and follow- 
up observations will be used to make a detailed investigation 
of the range and distribution of the intensity ratio for these 
two transitions. Combining this with our evolutionary stud- 
ies we will obtain unique insight into the changes in physical 
conditions responsible for the presence/absence of the dif- 
ferent maser methanol transitions. 

Our results relating to the evolution of star for- 
mation regions will be independently tested in the fu- 
ture through comparison between the unbiased sample of 
6.7-GHz methano l masers observed in the MMB survey 
l|Green et al.ll2009h . followup observations at 12.2-GHz and 
870 /xm dust sources observed as part of ATLASGAL 
(APEX Telescope Large Area Survey of the Galaxy j]]. 
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Predictor 


Deviance 


AIC 


LRT 


Pr(Chi) 


none 


455.47 


457.47 






Integrated 


429.24 


433.24 


26.23 


3.034e-07 


Peak 


417.15 


421.15 


38.32 


6.008e-10 


FWHM 


417.91 


421.91 


37.56 


8.881e-10 


Distance 


453.23 


457.23 


2.24 


0.1343038 


Mass 


445.21 


449.21 


10.26 


0.0013616 


Radius 


443.27 


447.27 


12.20 


0.0004782 


Density 


455.36 


459.36 


0.11 


0.7370057 



Table A2. Summary table for the Binomial regression model of 
1.2-mm dust clumps with and without 6.7-GHz methanol masers, 
showing for each predictor the estimated coefficient and the stan- 
dardised z-value and p-value for the test of the hypothesis that 
ft=0. 



Predictor 


Estimate 


Std. Error 


z value 


p-value 


Intercept 


-2.677e+00 


3.244e-01 


-8.251 


< 2e-16 


Integrated 


-2.952e-01 


7.338e-02 


-4.023 


5.75e-05 


Peak 


1.172e+00 


2.744e-01 


4.271 


1.94e-05 


FWHM 


2.994e-02 


6.427e-03 


4.659 


3.18e-06 


Density 


1.434e-06 


6.876e-07 


2.085 


0.0370 



APPENDIX A: DETAILS RELATING TO 
STATISTICAL ANALYSIS 

In this section we present the tables and box plots relating 
to the statistical analysis that is presented in Section 13.21 
All of the statistical analysis (including the production of 
box plots) w as completed using the R statistical analysis 
packag^f). See lBreen et al.l l|2007l . and references therein) for 
a more comprehensive description of the statistical methods 
adopted in our analysis. 

Box plots have been used as a way of presenting the data 
as they are a useful tool for displaying differences between 
populations. They can instantly hilight the distribution of 
the data as the separation between the different parts of the 
box have specific meanings. The solid line in the in each of 
these plots represents the median of the data. The box rep- 
resents the 25 th to the 75 th percentile, while the vertical line 
from the top of the box represents data from the 75 th per- 
centile to the maximum value and the vertical line from the 
bottom of the box represents data from the 25 th percentile 
to the minimum value. Extreme values, or outliers, are rep- 
resented by dots and are defined as points that exceed 1.5 
times the interquartile range in separation from either the 
25^ or 75th percentile. 



Table A3. Analysis of deviance table for the single term models, 
dust clumps with and without highly luminous 6.7-GHz methanol 
masers. 



property 


Deviance 


AIC 


LRT 


Pr(Chi) 


none 


140.006 


142.006 






Integrated 


139.777 


143.777 


0.229 


0.632381 


Peak 


139.690 


143.690 


0.316 


0.574197 


FWHM 


137.791 


141.791 


2.215 


0.136674 


Distance 


117.717 


121.717 


22.289 


2.345e-06 


Mass 


135.027 


139.027 


4.979 


0.025657 


Radius 


121.923 


125.923 


18.083 


2.115e-05 


Density 


131.441 


135.441 


8.564 


0.003428 


6.7 Maser flux 


95.731 


99.731 


44.275 


2.854e-ll 



Table A4. Summary table for the Binomial regression model 
of 1.2-mm dust clumps with luminous 6.7-GHz methanol masers 
and those with less luminous 6.7-GHz methanol masers, showing 
for each predictor the estimated coefficient and the standardised 
z-value and p-value for the test of the hypothesis that /3i=0. 



Predictor Estimate Std. Error z value p-value 

Intercept -3.183923 0.823230 -3.868 0.00011 

FWHM 0.017130 0.008611 1.989 0.04665 

Distance 0.488554 0.121787 4.012 6.03e-05 



2 R Development Core Team (2006). R: A language and en- 
vironment for statistical computing. R Foundation for Statis- 
tical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL 
htt p : / / www .R-project.org 
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Figure Al. Box plots of the 1.2-mm dust clump properties split 
in the categories of yes and no, referring to the presence or ab- 
sence of 6.7-GHz methanol maser emission respectively. There is 
a statistically significant difference between the two categories in 
the 1.2-mm dust clump integrated flux density, peak flux density, 
mass, radius and source FWHM. This information is the graphical 
display relating to the information in Table lATI 



Table A5. Analysis of deviance table for the single term models, 
dust clumps with 6.7-GHz methanol masers with and without 
12.2-GHz methanol. This analysis includes not only the 1.2-mm 
dust clump properties, but also both the peak flux density and 
luminosity of the associated 6.7-GHz methanol maser. 



Predictor 



Deviance 



AIC 



LRT Pr(Chi) 



none 


133.763 




135.763 




Integrated 


133.345 


137.345 


0.417 


0.51821 


Peak 


133.342 


137.342 


0.421 


0.51655 


FWHM 


131.321 


135.321 


2.442 


0.11814 


Distance 


133.588 


137.588 


0.175 


0.67568 


Mass 


133.559 


137.558 


0.205 


0.65090 


Radius 


133.573 


137.573 


1.190 


0.27529 


Density 


130.732 


134.732 


3.031 


0.08171 


6.7-GHz lum 


113.790 


117.790 


19.973 


7.853e-06 


6.7-GHz flux 


107.800 


111.800 


25.963 


3.481e-07 




Strong 6.7-GHz methanol maser 



yes 



Strong 6.7-GHz methanol maser 

Distance 



Strong 6.7-GHz methanol maser 



Mass 



no yes 
Strong 6.7-GHz methanol maser 

Radius 




Strong 6.7-GHz methanol maser 



H_2 Number Density 



Strong 6.7-GHz methanol maser 



no yes 
Strong 6.7-GHz methanol maser 

Figure A2. Box plots of the 1.2-mm dust clump properties split 
in the categories of yes and no, referring to the presence of 6.7- 
GHz methanol maser emission with high and low luminosities 
respectively. There is a statistically significant difference between 
the two categories in the 1.2-mm dust clump mass, radius and H2 
number density. This information is the graphical display relating 
to the information in Tabic [A3] 



Table A6. Summary table for the Binomial regression model of 
1.2-mm dust clumps with associated 6.7-GHz methanol masers 
with and without a 12.2-GHz methanol maser counterpart, show- 
ing for each predictor the estimated coefficient and the standard- 
ised z-value and p- value for the test of the hypothesis that 0i=O. 



Predictor 


Estimate 


Std. Error 


z value 


p- value 


Intercept 


-0.41058 


0.29388 


-1.397 


0.16238 


6.7-GHz flux 


0.03571 


0.01180 


3.027 


0.00247 
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Figure A3. Box plots of the 1.2-mm dust clump properties split 
in the categories of yes and no, referring to the presence or absence 
of 12.2-GHz methanol maser emission. There is a statistically 
significant difference between the two categories in the 1.2-mm 
dust clump H2 number density. This information is the graphical 
display relating to some of the information in Table 1X31 



Table A7. Analysis of deviance table for the single term mod- 
els, dust clumps with 6.7-GHz methanol masers with and with- 
out 12.2-GHz methanol. Data associated with 6.7-GHz with peak 
flux densities less than 9 Jy have been removed. The results of 
the single term models applied to the full data set are shown in 
Table lA5l 

Predictor Deviance AIC LRT Pr(Chi) 



none 


60.925 


62.925 






Integrated 


60.428 


64.428 


0.496 


0.481075 


Peak 


59.637 


63.637 


1.287 


0.256532 


FWHM 


59.856 


63.856 


1.069 


0.301177 


Distance 


59.343 


63.343 


1.582 


0.208483 


Mass 


60.924 


64.924 


0.0003 


0.986207 


Radius 


60.342 


64.342 


0.581 


0.445201 


Density 


55.748 


59.748 


5.177 


0.022894 


6.7-GHz Lum 


55.136 


59.136 


5.789 


0.016130 


6.7-GHz Flux 


53.484 


57.484 


7.440 


0.006377 




no yes 
12.2— GHz methanol maser presence 



6.7-GHz Luminosity 




no yes 
12.2— GHz methanol maser presence 

Figure A4. Box plots of the 6.7-GHz methanol maser properties 
split in the categories of yes and no, referring to the presence or 
absence of 12.2-GHz methanol maser emission. The differences 
in the two populations is highly significant in the case of both 
the 6.7-GHz methanol maser peak flux density and peak lumi- 
nosity. The three points with the most extreme values for either 
property have been removed to allow the distribution to be seen. 
This information is the graphical display relating to some of the 
information in Table lA"5l 



Table A8. Analysis of deviance table for the single term models, 
dust clumps with and without associated 8-GHz radio continuum. 



Predictor 


Deviance 


AIC 


LRT 


Pr(Chi) 


none 


365.19 


367.19 






Integrated 


340.56 


344.56 


24.63 


6.941e-07 


Peak 


330.03 


334.03 


35.17 


3.025e-09 


FWHM 


305.12 


309.12 


60.07 


9.145e-15 


Distance 


363.24 


367.24 


1.96 


0.161618 


Mass 


339.94 


343.94 


25.26 


5.023e-07 


Radius 


320.32 


324.32 


44.88 


2.096e-ll 


Density 


355.57 


359.57 


9.63 


0.001919 
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Figure A5. Box plots of the 1.2-mm dust clump properties split 
in the categories of yes and no, referring to the presence or absence 
of detectable radio continuum emission. There is a statistically 
significant difference between the two categories in all 1.2-mm 
dust clump properties with the exception of distance. The results 
of the single term models applied to the full data set are shown 
in Table |A8] 



